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ABSTRACT 

Using new UBVRIHa CCD photometric observations and the archival infrared and 
X-ray data, we have carried out a multi-wavelength study of a Perseus arm young 
galactic star cluster NGC 7419. An age of 22.5 ± 3.0 Myr and a distance of 3230t 33 o 
pc are derived for the cluster. Our photometric data indicates a higher value of color 
excess ratio E(U — B)/E(B — V) than the normal one. There is an evidence for 
mass segregation in this dynamically relaxed cluster and in the range 1.4 — 8.6M Q , 
the mass function slope is in agreement with the Salpeter value. Excess emissions in 
near-infrared and Ha support the existence of a young (< 2 Myr) stellar population 
of Herbig Ae/Be stars (> 3.0Mq) indicating a second episode of star formation in 
the cluster region. Using XMM-Newton observations, we found several X-ray sources 
in the cluster region but none of the Herbig Ae/Be stars is detected in X-rays. We 
compare the distribution of upper limits for Herbig Ae/Be stars with the X-ray dis- 
tribution functions of the T-Tauri and the Herbig Ae/Be stars from previous studies, 
and found that the X-ray emission level of these Herbig Ae/Be stars is not more 
than Lx ~ 5.2 x 10 30 erg s _1 , which is not significantly higher than for the T-Tauri 
stars. Therefore, X-ray emission from Herbig Ae/Be stars could be the result of ei- 
ther unresolved companion stars or a process similar to T-Tauri stars. We report an 
extended X-ray emission from the cluster region NGC 7419, with a total X-ray lumi- 
nosity estimate of ~ 1.8 x 10 31 erg s" 1 arcmin -2 . If the extended emission is due to 
unresolved emission from the point sources then we estimate ~288 T-Tauri stars in 
the cluster region each having X-ray luminosity ~ 1.0 x 10 30 erg s^ 1 . Investigation of 
dust attenuation and 12 CO emission map of a square degree region around the cluster 
indicates the presence of a foreground dust cloud which is most likely associated with 
the local arm star forming region (Sh2-154). This cloud harbors uniformly distributed 
pre-main-sequence stars (0.1 — 2.0M@), with no obvious trend of their distribution 
with either {H — K) excess or A v . This suggests that the star formation in this cloud 
depend mostly upon the primordial fragmentation. 

Key words: star clusters : young - star formation : primordial fragmentation - young 
stellar objects : Herbig Ae/Be, T Tauri, X-ray emission 



1 INTRODUCTION 

NGC 7419 (RAjaooo = 22 h 54 m 20 s , DEC.72000 = +60°48'54"; 
I = 109? 13, b — 1? 12), is a moderately populated young 
and heavily reddened galactic star cluster in Cepheus with 
a large number of Be stars. The cluster contains high mass 
(> IOMq); intermediate mass (2 — IOMq) and low-mass 
(< 2Mq) stars. It is therefore an ideal laboratory for the 
study of initial stellar mass distribution as well as duration 
of star formation process in a molecular cloud. Presence of 



statistically significant number of Herbig Ae/Be stars in the 
cluster makes it very attractive for understanding the forma- 
tion of these stars and origin of various atmospheric activi- 
ties like Ha emission and X-ray emission in them. However, 
to address these questions in detail, one would like to know 
accurate distance and age parameters of the cluster NGC 
7419, which is lacking despite a number of photometric and 
spectroscopic studies. This is mainly because of the fact that 
the cluster is heavily reddened in comparison to the nearby 
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clusters situated at the similar distances, and suffers from 
variable reddening. 

In order to determine cluster reddening reliably, accu- 
rate UBV broadband photometry of early type stars is es- 
sential. A comparison of the photometries available in the 
literature indicates that U B data, many have systematic cal- 
ibratio n error. For example, iBeauchamp. Moffat fc Drissenl 
1 19941 ) have mentioned that their color may have an offset 
of ~ 0.2 mag due to the calibration problems in U band. 
Their photometric observations have been carried out in the 
poor seeing (2'.'5 — 4'.'0) conditions. This will affect cluster 
photometric data particularly in the crowded regions. 

The main goals of present study are to determine the 
distance, age and its spread, and mass function (MF) of 
the cluster as accurate as possible. This will help us to un- 
derstand the star formation history of the cluster, and to 
investigate the X-ray emission properties of Herbig Ae/Be 
stars. Deep optical UBVRI observations (V ~ 22.0 mag), 
narrow band Ha photometric observations along with the 
Two Micron All Sky Survey (2MASS), Infrared All Sky Sur- 
vey (IRAS), Midcourse Space experiment (MSX) and XMM- 
Newton archival data are used to understand the X-ray emis- 
sion properties of Herbig Ae/Be stars, and the global sce- 
nario of star formation in the cluster NGC 7419 and its 
surr ounding regio n . 

iBlanco et al.l (1 19551 ) has reported the distance of 
this cluster as ~ 6 kpc based on the RI photomet- 
ric observations. A similar value for cluster d istance 
has also been obtained by Moffat fc VogtJ l| 19731 ) . How- 
ever, Ivan de Hulst et "ah I i|l954f ) has obtained a significantly 
smaller distance of 3.3 kpc. Using CCD data , a distance 
of 2.0 kpc and 2.3 kpc was estimated by Bhatt et al.l 
1 19931 ) and IBeauchamp. Moffat fc Drissenl ll 19941) . respec- 
tively . The age estimated bvlBhatt et all ll 19931 ) is ~ 40 Myr 
while IBeauchamp. Moffat fc Drissenl (Il994l ) have estimated 
a much younger age of ~ 14 Myr. Recent CCD observa- 
tions reported by ISubramaniam et alj |2006l) estimated its 
distance as 2.9 kpc and an age of 20-25 Myr. 

The paper describes optical observations and the deriva- 
tion of cluster parameters in [J2] and S}3] The near-infrared 
(NIR) data are dealt in S|4] while distribution of young stel- 
lar objects (YSOs), MF and mass segregation are given in 
Sj5] Sj6] and i}7] Finally, the X-ray data and its analysis (for 
the first time) are described in Sj8] followed by the summary 
and conclusions in SJU 



2 PHOTOMETRIC DATA 
2.1 Observations 

The optical observations of NGC 7419 were carried out using 
a thinned back- illuminated CCD camera mounted at f/13 
Cassegrain focus of the 104-cm Sampurnanand reflector tele- 
scope of Aryabhatta Research Institute of Observational Sci- 
ences, Nainital. A 24/im square pixel of the 2048 x 2048 size 
CCD detector corresponds to 0l.'38 and the entire chip cov- 
ers a field of about 13' x 13' on the sky. In order to improve 
the signal to noise ratio, observations were taken in binned 
mode of 2 x 2 pixel. The gain and readout noise of the CCD 
are 10 electron per analog-to-digital unit and 5.3 electron, 
respectively. The journal of optical observations is given 



Table 1. Journal of CCD obser vations of the c luster NGC 7419 
and the calibration region SA98 (Landolt 1992). 



Date(UT) 


Filter 


Exposure Time (s) 






(xno. of exposures) 


NGC 7419 


30 October 2005 


Ha 


300x2, 120x2, 1500x3 




Ha -cont 


300x2, 120x2, 900x3 


07 November 2005 


U 


300x3, 360x4 




B 


180x4, 240x4 




V 


180x4, 180x4 




R 


20x4, 120x4 




I 


20x4, 60x4 


08 November 2005 


U 


1800x1 




B 


1200x1 




V 


900x1 




R 


300x1 




I 


500x1 


25 October 2006 


U 


300x3 




B 


300x2 




V 


100x2 




R 


20x2 




I 


20x2 


SA 98 


25 October 2006 


U 


360x11 




B 


150x10 




V 


60x10 




R 


25x10 




1 


25x10 



in Table \T\ Broad band Johnson UBV, Cousins RI and, 
narrow band Ha line (A = 656.5 nm) and Ha continuum 
(A = 665 nm) filters were used for observations. The nar- 
row band filters had a full width at half maximum of 8 nm. 
Several bias and twilight flat field frames in all the filters 
were taken to clean the images. Multiple long and short ex- 
posures w ere obtained fo r the cluster region. We observed 
12 stars in iLandoltl l|l992l ) standards field (SA98) covering a 
range in brightness (11.93 < V < 15.90) as well as in color 
(0.157 < (B -V) < 2.192) for calibrating the cluster ob- 
ser vations. The ai rmass range is covered from 1.1 to 2.0 for 
the lLandoltl (|l992l ) standards, which was used for extinction 
determinations. 

The photometric CCD data were reduced using the 
IRAFQ and ESO MIDAS^ data reduction packages. Pho- 
tometry of the bias-subtracted and flat-fielded CC D frames 
was carrie d out using DAOPHOT-II software l|Stetsonl 
Il987l , ll992l) . Magnitude of the stars obtained from different 
frames were averaged separately for short and long expo- 
sures. When brighter stars were saturated on deep exposure 
frames, their magnitudes were taken only from unsaturated 
short exposure frames. We used DAOGROW program for 



1 IRAF - Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science 
Foundation ( http: / /iraf.noao.edu ). 

2 MIDAS — Munich Image Data Analysis System is developed 
and maintained by ESO, the European Southern Observatory. 
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Table 2 . The mean photometric errors (<r) in magnitude are given 
as a function of brightness range in the cluster region. 



Magnitude 


<J,r 




<n- 




07 


range 












< 14 


0.009 


0.002 


0.003 


0.005 


0.007 


14 - 15 


0.007 


0.005 


0.003 


0.005 


0.004 


15 - 16 


0.021 


0.010 


0.013 


0.012 


0.025 


16 - 17 


0.030 


0.015 


0.010 


0.011 


0.023 


17- 18 


0.069 


0.018 


0.010 


0.012 


0.026 


18 - 19 


0.131 


0.029 


0.014 


0.011 


0.031 


19 - 20 




0.070 


0.017 


0.013 


0.026 


20 - 21 




0.173 


0.031 


0.019 


0.031 


21 - 22 






0.082 


0.036 


0.034 


22 - 23 






0.196 


0.061 


0.047 



0.2 

AU o 



l_ ( 


*) SA 98 
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Figure 1. Difference in UBVRI magnitude measurements as a 
function of V magnitude are shown. T he left panels (a) show 
the residuals of standard mag nitude from lLandoltl dl992f ) and the 
present transformed magnitudes, while the middle (b) and right 
(c) panels show magnitude difference of secondary standards as 
generated on 25 October 2006 to that of on 07 November 2005 
and 08 November 2005, respectively. 



construction of an aperture growth curve required for de- 
termining the difference between aperture and profile fitting 
magnitudes. 



2.2 Calibrations 

The instrumental magnitudes were converted into the stan- 
dard system using le a st squ are linear regression procedure 
outlined by IStetsonl (| 1992T ) . The photometric calibration 
equations used are: 

u = f/+(0.588±0.015)X+(--0.002±0.004)((7-B) + (6.933±0.006) 
b = B+(0.347±0.013)X+(-0.035±0.005)(B- V) + (4.771±0.007) 
v = V+(0.159±0.012)X+(-0.041±0.004)(V-i") + (4.320±0.007) 
r = iJ+(0.139±0.009)X+(-0.056±0.007)(V-i?)+(4.234±0.006) 
i = 7+(0.104±0.007)X + (-0.048±0.003)(V-/) + (4.767±0.005) 



where U, B, V, R and / are the standard magnitudes; u, b, v, r 
and i are the instrumental magnitudes obtained after time 
and aperture correction; and X is the airmass. We have ig- 
nored the second-order color correction terms as they are 
generally small in comparison to the internal photometric 
errors as given by DAOPHOT. The photometric errors as 
a function of brightness range are given in Table [2] It can 
be seen that the errors become large (> O.lmag) for stars 
fainter than V « 22 mag, and hence the measurements be- 
yond this magnitude are less reliable and are not considered 
in further analysis. At V band, we could detect 1817 stars in 
13' x 13' region and their photometric magnitudes are given 
in Table [3] (available only in electronic form). The difference 
of calibrated mag nitudes derived using above transforma- 
tion to that of the lLandoltl (|l992h magnitudes are plotted in 
Figure [1] We generated secondary standards in the cluster 
field using data of 25th October 2006 to standardise the data 
observed on 7th and 8th of November 2005. We, therefore 
also plot differences in magnitudes of the secondary stan- 
dards on these dates. No systematic effect has been seen 
in the residuals from night-to-night and its distribution is 
random in nature with a typical accuracy of ~ 0.03 mag in 
UBVRI band. 

Figure [2] shows a comparison of the present CCD pho- 
tometry wit h the p r eviously reported CCD photom e try by 
iBhatt et all dl993T), IBeauchamp, Moffat fc Drissenl (|l994T ) 
and ISubramaniam et al.l (I2006T ). The difference A fpresent- 
literature) is plotted as a function of V magnitude and a 
detailed statistics are given in Table [4] (available only in 
electronic form). In comparison to th e photometric data of 
IBeauchamp. Moffat fc Drissenl (|l994T l. the (B - V) colors 
obtained by us are bluer by ~ 0.25 mag, (U — B) col- 
ors are redder by ~ 0.25 mag and the V magnitude is 
faint er by ~ 0.05 mag. Our p hoto metry is in a greem ent 
with lSubramaniam etall l^OOfj ) and IBhatt et all (|1993T l. It 
is worth pointing out that our data is ~ 2.0 mag de eper from 
ISubramaniam et al.l (|2006t ) and lBhatt et all l| 1993ft . i.e., V ~ 
22.0 mag. 



2.3 Completeness of the data 

The completeness of the data used in the present work was 
estimated using the ADDSTAR routine of the DAOPHOT 
II. In brief, the method consists of randomly adding ar- 
tificial stars (about 10-15% of the originally detected 
stars) of known magnitudes and positions into the orig- 
inal frames. The frames were re-reduced using the same 
procedure used for the original frames. The ratio of 
the number of recovered stars to those added in each 
magnitude interval gives the completeness factor (CF) 
as a function of magnitude. The CF was obtained us- 
ing the stars which were recovered in both V and I 
pass-bands . The detailed procedures have been outlined 
elsewhere (ISapr fc Richtlerl Il99ll; ISagar fc Griffiths] Il998l ; 



iPandev et al.N200ll ; iNilakshi fc Sagarll2002l ). The CF 
function of V magnitude is given in Table [S] 



2.4 Ha photometric data 

To identify emission line stars, we use a-index (rri6565 — 
m-6650) parameter, where TO6565 and meeso are the magni- 
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Table 3. UBVRI and Ha photometric data of the sample stars in 13'xl3' region. X and Y positions of stars in the CCD are converted 
into RAj2ooo an d DECj2000 using the Guide Star Catalogue II (CSC 2.2, 2001). The complete table is available in electronic form only. 



ID 


RA .J2000 
(h m s) 


DEC J2Q00 


X 

(pixel) 


Y 

(pixel) 


(" 


u 

iag 


) 




C* 


B 

iag) 




(* 


V 

lag) 








R 


) 






(« 


I 

iag 


) 


Ha 


na 


dcx 

I) 


1 


22 54 06.38 


+ 60 52 13.30 


818.710 


640.714 


12.999 


± 


0.002 


12 


481 


± 0.002 


11 


707 


± 


009 


11 


321 


± 





012 


11 


003 


± 


0.002 


0.032 


± 


0.008 


2 


22 54 02.22 


+ 60 54 57.20 


1032.410 


691.088 


13.518 


= 


0.002 


12 


964 


± 0.002 


12 


156 


± 


009 


11 


761 


= 





011 


11 


381 


= 


0.004 


0.038 


± 


0.006 


3 


22 53 30.17 


+ 60 53 38.60 


914.148 


994.003 


13.964 


= 


0.003 


13 


492 


± 0.002 


12 


740 


± 


002 


12 


321 


= 





010 


11 


946 


= 


0.004 


0.022 


± 


0.020 

























































Tabl e 4. A comparison of present U BVRI CCD photometry with that present in the literature - viz Bha tt et al 
Beauchamp, Moffat & Drisscn (1994) and Subramaniam et al. (2006). The complete table is available in electronic form only. 



V range < AV > 


< A(B -V)> 


< A(U-B) > 


< A(V -R)> 


< A{V-I) > 


Mean ± a 


Mean ± a 


Mean ± <r 


Mean ± a 


Mean ± a 


Bhatt et al. (4993 ) 










< 14.0 0.025 ± 0.030 (8) 


0.060 ± 0.070 (6) 


0.030 ± 0.083 (4) 


-0.022 ± 0.065 (7) 


-0.044 ± 0.034 (6) 














ihamp at al. (1994) 



•- 'i'*^«»:>^ 



14 16 18 20 14 16 18 20 14 16 18 20 

V 

Figure 2. Comparison of present CCD photometry with that 
available in the literature. The difference A (present-literature) 
is in magnitude. The dashed line drawn in each panel represents 
zero difference. 



Table 5. Completeness Factor of the photometric data as a func- 
tion of brightness and location. 



Magnitude range 




NGC 7419 


Field region 


(V mag) 


r < 1 


1' < r < 3.5' 


3.'5 < r < 5' 


13-14 


1.00 


1.00 


1.00 


14-15 


1.00 


1.00 


1.00 


15-16 


1.00 


1.00 


1.00 


16-17 


1.00 


0.98 


0.99 


17-18 


0.89 


0.97 


0.98 


18-19 


0.57 


0.93 


0.96 


19-20 


0.59 


0.89 


0.97 


20-21 


0.60 


0.86 


0.93 


21-22 


0.15 


0.38 


0.47 



iii 

~ Emission Line Stars 

1 

i 


\ 1 ' : 

[piiJ : 


— * I . * " 




Absorption— Line Stars " • 


0.8 1.2 


1.6 

(R I) 



-0.8 
-0.7 
-0.6 
-0.5 
-0.4 
-0.3 
-0.2 
-0.1 
0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 



Figure 3. The Q-index as a function of color (R — I). The dashed 
line represents the 3cr limits from a mean o-index of 0.096±0.116 
mag. Emission-line stars having o-index < —0.252 mag are shown 
with error bars. The absorption line stars denote a-index > 0.444 
mag. 



tude of a star in the Ha and Ha -continuum filter, respec- 
tively. In an area of 13' x 13', we could measure a-index for 
1065 stars and the same is plotted against (R — I) colors in 
Figure [3] We estimated a mean a-index iteratively for these 
stars by 3-sigma clipping. Seven such iterations resulted in 
a constant mean value of 0.097 mag and a RMS scatter 
(a) of 0.116 mag. Stars having a-index within 3<r are there- 
fore considered as having normal strength at Ha feature. 
A star is considered as Ha emitter if a < —0.25 mag and 
Ha absorber if a > 0.44 mag. The observed characteristics 
of stars with Ha emission are given in Table [6] Ha emission 
stars reported only by others in the literature are also listed 
in Table [75] This could be due to the variable nature of 
Ha emission from the stars. 

We detect 29 Ha emission-line stars, and along 
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Table 6. Photometric properties of 44 Ha emission-line stars in 13' X 13' region around NGC 7419 are tabulated. Columns 1-9 are 
taken from the present work. Column (10) carry comments on the stars identified in the present and pr evious wo rk - 'Yes' stands for 
detect ion and 'N o' for non-detection in the present work, 'BMD' forlBcauchamp, Moffat & Drissen l|l994h . 'PK' for lPigulski fc Kopackil 
(2000), ' Sub' for lSub ramania m et al. I j2006l) along with their identification in braces. Stars located outside the cluster region (see i|2.4H 
are indicated with 'field'. 



ID RAj2ooo DECj2ooo X Y 

(h m s) (° ' ") (pixel) (pixel) 
(1) (2) (3) (4) (5) 



V R a-index R — I Remark 

(mag) (mag) (mag) (mag) 
(6) (7) (8) (9) (10) 



12 


22 


54 


1 1 


60 


18 


39 


541 


2 


549.5 


13 


75 


12 


63 


-0 


36 


+ 1 


08 


Yes, BMD(389), PK, Sub(M) 


15 


22 


54 


39 


60 


47 


24 


454 


5 


310.5 


13 


93 


12 


96 


-0 


01 


+0 


97 


No, BMD(1129), PK 


18 


22 


54 


26 


60 


49 


02 


576 





442.9 


14 


12 


13 


15 


-0 


10 


+0 


96 


No, BMD(781), PK, Sub(C) 


19 


22 


54 


20 


60 


49 


52 


639 


5 


496.0 


14 


14 


13 


05 


-0 


10 


+1 


05 


Yes, BMD (620), PK, Sub(G) 


20 


22 


54 


15 


60 


19 


50 


634 


7 


545.2 


14 


11 


13 


11 


-0 


18 


+1 


03 


No, BMD(417), PK, Sub(J) 


21 


22 


54 


2.3 


60 


•50 


04 


656 


I 


473.7 


14 


20 


13 


07 


-0 


38 


+1 


10 


Yes, BMD(702), PK, Sub(4) 


23 


22 


54 


2.3 


60 


49 


53 


642 


3 


475.6 


14 


23 


13 


09 


-0 


56 


+1 


07 


Yes, BMD(692), PK, Sub(3) 


27 


22 


54 


20 


60 


48 


36 


539 





499.1 


14 


36 


13 


33 


-0 


21 


+1 


02 


No, BMD(585), PK, Sub(I) 


29 


22 


54 


37 


60 


48 


36 


547 


6 


335.7 


14 


39 


13 


23 


-0 


50 


+1 


10 


Yes, BMD(1076), PK, Sub(A) 


38 


22 


54 


27 


60 


48 


52 


564 





428.1 


14 


88 


13 


76 


-0 


54 


+1 


04 


Yes, BMD(831), PK, Sub(B) 


39 


22 


54 


20 


60 


48 


54 


563 


3 


492.1 


14 


93 


13 


92 


-0 


05 


+1 


04 


No, BMD(621), PK, Sub(K) 


18 


22 


54 


18 


60 


48 


57 


566 


2 


517.2 


15 


14 


14 


08 


-0 


50 


+1 


01 


Yes, BMD(504), PK, Sub(L) 


57 


22 


54 


26 


60 


19 


05 


581 


2 


435.0 


15 


37 


11 


30 


-0 


24 


+0 


81 


No, BMD(815), PK, Sub(2) 


59 


22 


54 


08 


60 


50 


23 


673 


6 


621.9 


15 


44 


11 


18 


-0 


55 


+1 


20 


Yes, BMD(239), PK, Sub(5) 


64 


22 


54 


16 


60 


47 


49 


476 





532.4 


15 


54 


14 


32 


-0 


57 


+1 


16 


Yes, BMD(427), PK, Sub(N) 


66 


22 


54 


14 


60 


46 


20 


358 





543.2 


15 


56 


14 


51 


-0 


22 


+1 


13 


No, BMD(375), Sub(P) 


69 


22 


54 


24 


60 


49 


31 


613 


6 


463.4 


15 


70 


14 


69 


-0 


11 


+0 


97 


Yes, BMD(728), PK, Sub(D) 


71 


22 


54 


20 


60 


48 


51 


559 


3 


500.3 


15 


71 


11 


76 


-0 


08 


+0 


96 


No, BMD(582), PK, Sub(H) 


73 


22 


54 


32 


60 


47 


57 


494 


2 


374.4 


15 


72 


14 


65 


+0 


04 


+1 


15 


No, BMD(967), PK 


74 


22 


54 


10 


60 


49 


39 


617 


3 


595.2 


15 


73 


14 


67 


+0 


06 


+0 


59 


No, BMD(290), PK 


75 


22 


54 


29 


60 


19 


08 


586 


1 


408.3 


15 


78 


14 


68 


-0 


52 


+1 


06 


Yes, BMD(884), PK, Sub(l) 


81 


22 


54 


17 


60 


18 


23 


520 


3 


526.5 


15 


95 


14 


98 


-0 


02 


+0 


99 


No, BMD(451), PK 


85 


22 


54 


17 


60 


48 


18 


514 


4 


524.8 


15 


97 


14 


93 


+0 


05 


+1 


06 


No, BMD(458), PK 


89 


22 


54 


15 


60 


•51 


22 


755 





556.0 


16 


01 


14 


92 


_0 


55 


-1-1 


06 


Yes, Sub(Q) 


91 


22 


54 


13 


60 


50 


09 


658 





572.1 


16 


06 


15 


00 


-0 


01 


+1 


13 


No, BMD(351), PK 


95 


22 


54 


17 


60 


51 


37 


776 


5 


535.6 


16 


11 


15 


10 


-0 


17 


+1 


01 


Yes, Sub(R) 


110 


22 


54 


26 


60 


47 


57 


491 


7 


435.1 


16 


35 


15 


28 


-0 


36 


+1 


04 


Yes, BMD(795), PK, Sub(6), 


126 


22 


54 


12 


60 


48 


33 


532 


6 


569.7 


16 


50 


15 


43 


-0 


01 


+1 


13 


No, BMD (340), PK 


127 


22 


54 


19 


60 


48 


14 


509 


7 


506.9 


16 


50 


15 


45 


-0 


08 


+1 


11 


No, BMD(551), PK 


139 


22 


54 


21 


60 


47 


36 


462 


8 


450.1 


16 


65 


15 


43 


-0 


51 


+1 


15 


Yes, BMD(745), Sub(E) 


115 


22 


54 


07 


60 


48 


17 


508 


8 


618.9 


16 


71 


15 


33 


-0 


29 


+1 


33 


Yes 


181 


22 


54 


51 


60 


48 


27 


543 


4 


199.9 


17 


07 


15 


95 


-0 


28 


+1 


21 


Yes, field 


218 


22 


54 


24 


60 


47 


01 


416 


6 


447.5 


17 


32 


16 


17 


-0 


40 


+1 


12 


Yes, BMD(741), PK, Sub(F) 


321 


22 


54 


1.3 


60 


50 


27 


681 


9 


572.7 


18 


13 


16 


92 


-0 


32 


+1 


23 


Yes, BMD (352), 


438 


22 


54 


38 


60 


52 


10 


829 


5 


331.6 


18 


81 


17 


55 


-0 


30 


+1 


36 


Yes, field 


_t 


22 


53 


3.3 


60 


45 


51 


300 


2 


938.4 


18 


87 


17 


27 


-0 


28 


+2 


11 


Yes, field 


655 


22 


54 


56 


60 


47 


41 


485 


2 


148.1 


19 


53 


18 


14 


-0 


16 


+1 


58 


Yes, field 


703 


22 


53 


37 


60 


53 


22 


895 


1 


925.5 


19 


65 


18 


45 


-0 


26 


+1 


24 


Yes, field 


714 


22 


54 


10 


60 


49 


24 


612 


6 


306.4 


19 


68 


18 


40 


-0 


35 


+1 


11 


Yes, BMD(1157) 


716 


22 


54 


57 


60 


53 


38 


954 


5 


157.8 


19 


68 


18 


37 


-0 


32 


+1 


42 


Yes, field 


791 


22 


55 


05 


60 


53 


02 


911 





82.3 


19 


90 


18 


52 


-0 


13 


+1 


51 


Yes, field 


801 


22 


54 


34 


60 


51 


19 


760 


3 


366.9 


19 


92 


18 


59 


-0 


32 


+1 


47 


Yes 


808 


22 


54 


10 


60 


18 


31 


542 


9 


303.1 


19 


91 


18 


70 


-0 


26 


+1 


26 


Yes, BMD(1153) 


831 


22 


54 


13 


60 


45 


59 


329 


3 


553.2 


20 


01 


18 


45 


-0 


18 


+1 


78 


Yes 



Notes: * — Error in Vmagnitude is greater than 0.1 mag. 



with 15 more reported by IPigulski fc Kopackil {2000) 
and ISubramaniam et al.l l)200fj| ). the total number of 
Ha emission-line stars in the 13' x 13' region around the 
cluster have now increased to 44. Of these, thirteen are 
newly identified in the present work. Though some late-type 
dwarfs do show strong chromospheric Ha in emission, but a 
few false detection may not be ruled out. As a non-emission 
late-type star generally has a series of strong metal oxide ab- 
sorption lines, such as TiO in its spectrum and this may lead 



to an under-estimation o f their continuum, hence mak ing the 
star an Ha emitter (e.g. iHuang. Chen fc H siao 200l). 



3 PROPERTIES OF THE CLUSTER 
3.1 Radial density profile 

We estimate the cluster center iteratively by calculating 
average X and Y positions of stars with V < 18.0 mag 
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2 f 4 r = 0.92±0. 14 arcmin 

f b = 0.41±0.07 stars arcmin" : 



(V< 17 mag) 




r c = 1.09±0.03 arcmin 
2 ^ f b = 2.87±0.19 stars arcmin" 5 



(V< 20 mag) 




r e = O.B4±0.10 arcmin 
1.5 h f b = 3.17±0.18 stars arcmin-' 



(2MASS) 




2 4 
Radius (arcmin) 

Figure 4. Projected radial stellar density profile of NGC 7419. 
Dashed lines represent 3<r levels above the field sta r dens i ty an d 
solid curve shows best fit to the empirical model of iKina (119620 . 



within 80 pixels from eye estimated center, until it converged 
to a constant value. The coordinates of the cluster center 
are found to be RA J20 oo = 22 h 54 m 18 s and DEC.72000 = 
+60°48'58". Using this method, a typical error expected in 
locating the center is 5". To estimate extent of the cluster, 
we divide the cluster region into a number of concentric cir- 
cles with an annulus of width 30'.'5. The projected radial 
stellar density profile is plotted in Figure [4] for three magni- 
tude limits - e.g. V—17, 18, 20 mag. The error bars denote 
Poisson statistics. The p( r) is parameterized as (|King|| 19621 ; 
iKaluznv fc Udalski|[l992) 



p(r) oc 



fa 



(1) 



1 + (r/r c ) 2 

where r c and /o are the core radius of the cluster and central 
star density, respectively. As the observed area (13' x 13') is 
large compared to the cluster diameter, we estimate /b us- 
ing outer region (r > 5') of the cluster and to derive /o and 
r c , we p erformed a Leven berg-Marquardt non-linear fitting 
routine IjPress et al.ll 1992T ) to equation [1] The fitted param- 
eters are shown in Figure [4] We define extent of the cluster 
(r c i) as the radius at which p(r) equals 3<j level above the 
field star density. Using optical data, we have estimated the 
core radius of about 1' and cluster radius of 4.'0±0!5 for the 
cluster. This is also supported by the radial density profile 
fits at 2MASS K-ba nd data. The catalogue of open cluster 
by iDias et ail l|2002i ) have been reported the radius of the 
cluster as 2'.5. To reduce the effect of field star contamina- 
tion, we therefore, consider lower limit of radius as 3^5 for 
further analysis. 



3.2 Color-magnitude diagrams 

The V, (U-B); V, (B — V) and V, {V-I) color-magnitude di- 
agrams (CMDs) of the cluster region along with V, (V — I) 
CMD of the field region are shown in Figure [5] Morphol- 
ogy of the CMDs is typical of a young-age open star clus- 



ter and the main sequence extends down to V ~ 22.0 
mag in V, (V — J) CMD. It is seen that the field star con- 
tamination becomes more evident for stars with V > 18 
mag, therefore, a statistical criteria was used to remove the 
field star contamination, follow ing procedure adopted by 
ISandhu. Pandev fc Sagarl (|2003l ). For this purpose, we as- 
sume a uniform distribution of field stars in an area outside 
the cluster region (r > 3!5). The luminosity function (LF) 
derived from the V, (V — I) CMD of the field region was 
subtracted from the LF of the cluster region. For a ran- 
domly selected star in the V, (V — I) CMD of field region, 
the nearest star in the V, (V — I) CMD of cluster region 
within V ± 0.25 mag and (V-I) ± 0.125 mag of the 
field star was removed. The necessary corrections for CF 
have also been considered in statistical subtraction. Statis- 
tical subtraction of data have been done till V w 21.0 mag 
only as they have CF value > 0.5. The V, (V — I) diagram 
of statistically cleaned cluster sample is shown in Figure [S] 
Further, to select main-sequence members more reliably, we 
define blue and red envelope in statistically cleaned CMD 
and the same is shown in the left panel of Figure [U] 



3.3 Interstellar reddening 

3. 3. 1 Law of interstellar reddening 

We investigated the nature of interstellar reddening towards 
the cluster d irection using th e color excess ratio method as 
described by I Johnson! (|l968T ). We select stars with spectral 
types earlier than AO by applying criterion of V < 17.7 and 
1.35 < (B — V) < 1.70 on probable members of the clus- 
ter ( H3.2p . From this list we exclude stars showing either 
Ha emission (see £|2.4[) or NIR excess (seeQ as their red- 
deni ng properties are lik ely to be different from the normal 
stars Kumar et al.l (120041). We estimate i ntrinsic colors using 
Q-method fcf. Ijohnson fc Morgan! Il953l) and iteratively es- 
timate reddening free parameter Q [= (U — B) — X(B — V), 
where X = E(U - B)/E(B - V)], till the color excess ra- 
tios become constant within the photometric errors. For the 
first iteration, we assume X — 0.72. The color excesses 
are determined using intrinsic colors as derived from the 
MKK spectral t ype-l uminosity class color relation given by 
ICaldwell et all Il993l) for ( U - B ), (B - V), (V - R) and 
{V - I) and bv lKoornneefl i| 19831 ) for (V — J), (V — H) and 
(V — K). Mean values of color excess ratios derived for the 
cluster are given in Table [7] and except X, all other color 
excess ratios are in agre ement with that expected for nor - 
mal interstellar matter l|Cardelli. Clayton fc Mathislll989T ). 
The value of X (~ 0.90 ± 0.08) is significantly greater than 
the normal value of ~ 0.72, which implies that the aver- 
age dust grain sizes are smaller than average. The ratio 
of total-to-selective extinction is esti mated using the rela- 
tion R v = 1.1E(V - K)/E(B - V) (|Whittet fc van Bredal 
Il98d ) and it is found to be 3.2 ± 0.1. 

We have also determined the value of color excess ratios 
using our UBVRI photometry and spectral cla ssification of 
six stars (see TableEl) available in the literat ure (|Caron et al.l 
20031 ; iBeauchamp. Moffat fc Drissenl 1 19941 ) . Intrinsic colors 
3 are read from Iwegnerl <|l994h for (V - R) and (V - I), 



The intrinsic colors of Johnson RI in IWegnerl lll994) are con- 
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0123 1 2 3 1234 1234 1234 

(U-B) (B-V) (VI) (V-I) (V-I) 

Figure 5. Color-magnitude diagrams for cluster regions are shown in panels (a) V, (V — B) (b) V, (B — V) and (c) V, (V — I). Panel (d) 
is for field region and the statistically cleaned V, (V — I) CMD for cluster is shown in panel (e). Theoretical isochrones from GRD02 are 
shown with continuous lines for log(age)= 6.6, 7.3 and 7.4 Myr in the panels (a), (b), (c) and (c). 



Table 7. Observed color exces s ratios in the direction of NGC 7419 are listed along with the excess ratios valid for normal interstellr 
reddening law, R v =3.l llCardelli. Clayton k. Mathislll989T) . 



Object 


E(U-B) 
E(B-V) 


E(V-R) 
E(B-V) 


E{V-I) 
E(B-V) 


E(U-B) 
E(V-J) 


E(B-V) 
E(V-J) 


E(V-R) 
E(V-J) 


E(V-I) 
E(V-J) 


E(V-H) 
E(V-J) 


E{V-K) 
E(V-J) 


NGC 7419 


0.90±0.09 


0.62±0.02 


1.29±0.04 


0.37±0.04 


0.41±0.01 


0.25±0.01 


0.52±0.01 


1.14±0.01 


1.20±0.02 


Normal value 


0.72 


0.62 


1.25 


0.32 


0.43 


0.27 


0.56 


1.13 


1.21 



from iFitzgeraldl (|l970l ') for (U-B). The mean value of X 
(= 0.88 ±0.06) using this method is in good agreement with 
the one derived from iterative Q-method. Further, the excess 
ratios at other wavebands, and mean value of Rv (= 3.16 ± 
0.10) are also close to the one derived using Q-method (see 
Table 0. 

Higher value of color excess ratio X, as found above, 
is also supported by the V, (U - B) C MD (Figure 0, in 
which the isochrones l|Girardi et al.ll20o3 . hereafter GRD02) 
are found to fit well with X = 0.90. Such anomalous val- 
ues of X are not uncommon for the interstellar matter in 
the Milky way and a higher value of X(> 0.72) has also 
been reported along lines of sight of ot her star forming clus - 
ters, for example, NGC 8 69 (X = 0.90;|Pandev et all 2003h ; 
Markarian 50 (X = 0.81; lBaume. Vazquez fc Carrardl2004t) 
and NGC 7510 (X = 0.78; iBarbon fc Hassanlll996l ) How- 
ever, they report normal value for other color excesses in 
these clusters. For further analysis, we have therefore used 
the normal value of color excess ratios except for the value 
E(U - B)/E(B - V), which is adopted as 0.90. 

3.3.2 Determination of reddening 

Figure [7] shows observed color-color, (U — B),(B — V), 
diagram for the cluster region. To determine reddening 
E(B — V), the in trinsic colors for main -sequence stars with 
solar metallicity l|Schmidt-Kalerl [l982f ) are reddened along 





(V i) 



(V-I) 



Figure 6. Statistically cleaned V, (V - I) CMD for NGC 7419. 
Left panel: Solid lines confine the blueward and redward bound- 
ary of the main-sequence. Emission line stars are encircled. Right 
panel: Solid curves show the be st fit Padova isochrone obtained 
using T 2 -minimization method l|Navlor fc JeffriedExjh for fixed 
z=0.019 and age =25 Myr with distance as free parameter. Hor- 
izontal bars denote errors in (V — I). See 93. 51 for further details. 



verted into Cousins RI using iBesselll (Il979h transformation rela- 
tions. 
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Table 8. Giant stars with known spectral types in the literature (jCaron ct al. 2003; Bca uchamp. Moffat fc DrissedH994h . Color excess 
ratios are derived using present data. 



ID 


BMD ID 


Spectral Type 


E(B-V) 


E(U-B) 
E(B-V) 


E(V~R) 
E(B-V) 


E(V-I) 
E(B-V) 


E(V-K) 
E(B-V) 




K) 

-V) 


8 


687 


B2.5 II 


1.741 


0.883 


0.641 


1.283 


2.924 


3.216 




13 


473 


B2.5 III 


1.648 


0.891 


0.623 


1.245 


2.754 


3.030 




15 


1129 


B4.0 III 


1.669 


0.752 


0.639 


1.266 


2.833 


3.116 




1(5 


190 


BO.O III 


1.999 


0.924 


0.642 


1.267 


2.837 


3.120 




20 


417 


Bl/2 III 


1.845 


0.912 


0.638 


1.263 


2.903 


3.193 




12 


389 


Bl/2 III 


1.789 


0.896 


0.670 


1.313 


3.014 


3.315 






(B V) 

Figure 7. The (B — V), (U — B) color-color diagram for stars in 
the cluster region. Open circles represent stars with Ha emission. 
Solid lines are intrinsic main-sequence reddened along the red- 
dening line with E(B — V) = 1.5 mag and 1.9 mag, while the 
dotted line is for 1.7 mag. 



the reddening vector (X — 0.9). This reddened main- 
sequence suggests that E(B—V) for the cluster region varies 
from 1.5 to 1.9 mag, with a visual mean value of about 
1.7 mag. This reddening variation (~ 0.4 mag) for cluster 
members is well above the typical width of intrinsic main- 
sequence stars i.e. < 0.11 mag, which is caused usually due 
to photometric ac curacy and the presence of binary stars 
l|Sagar et all I1987T ). It therefore indicates the presence of 
variable red dening within the cluster r egion, which is also 
reported by ISubramaniam et al.l (|2Q0d) . Few emission-line 
stars (open circles) are found to have E(B — V) > 1.9 mag 
and these may either be embedded deep in the parent cloud 
or be surrounded with circumstellar matter. Using iterative 
Q-method f £|3.3.ip , we derive mean value of E(B — V) as 
1.7 ±0.2 mag, which is same as estimated visually by color- 
color diagram, and hence we adopt this value in our further 
analysis. 



3.4 Turn-off age of the cluster 

Turn-off age of the cluster is determined by comparing 
the theoretical stellar evolutionary models, the Padova 



isochrones GRD02 for solar metallicity Z=0.019 with ob- 
served CMDs (see Figure [5}. Theoretical isochrones for 
ages 4, 20 and 25 Myr, corrected for the mean reddening 
E(B — V) — 1.7 mag, are shown with solid lines and are visu- 
ally fitted to the bluest envelope of CMDs consisting of prob- 
able members as selected in £13.21 It is found that X = 0.9 
results in a better isochrone fit in V, (U — B) CMD than 
the normal value, X = 0.72 f i|3.3.1| l. Locations of the red 



supergiants (jCaron et al.ll2003h BMD 921 (M2.5Iab), BMD 
696 (M2.5Iab) and BMD 435 (M2.5Iab) are found to be 
consistent with the isochrones, while BMD 139 (M3.5I) and 
BMD 950 (M7.5), both pulsating supergiants were prone to 
produce more scatter in distance and age estimate of NGC 
7419. Thus, using the visual-fit method, the present data 
finds the turn-off age of the cluster to lie between 19 Myr to 
25 Myr and hence we assign an age of 22.5 ± 2.5 Myr. 

We also estimate age using morphological age metho d 
l|Phelps. Janes fc Montgometrvl [19941 : iPandev et ail ll997T ). 
which employes color-index parameter (CIP), defined as the 
difference in the color index between the blue turnoff point 
of the main-sequence and the color at the base of the red 
giant branch. For NGC 7419, CIP is calculated as 1.654 mag 
and it represents the age between 20 and 25 Myrs, which is 
in good agreement with our previous va lue. The present age 
estima te is similar to the one derived bvlSubramaniam et al.l 
( 2006) less than the 50 Myr given bvlBhatt et alj (ll993Tl and 
older t han 14 Myr derived bv lBeauchamp. Moffat fc Drissen] 
(|l994h . 



3.5 Distance to the cluster 

The zero-age main-sequence fitting procedure (Figure [5} 
is used to derive the distance of the cluster and the dis- 
tance modulus is estimated as 12.45 ± 0.20 mag which cor- 
responds to a distance of 3100 ± 290 pc. The distance de- 
termination is further performed using a maximum likeli- 
hood r 2 -minimization method which employes fitting two- 
dimensional d istributions to stellar d ata in color mag- 
nitude space dNavlor fc Jeffries! I2006J : iMavne et all 120071 ; 
iJeffries et al. I l2007h . This takes into account the effects of 
binary population as well as observed photometric errors. 
Assuming an age of ~ 22.5 Myr and a binary fraction of 0.5, 
the best-fit Padova isochrone with solar metallicity yields a 
distance modulus of 12.55 mag (Figure ltTl H3.2|l , Emission- 
line stars were excluded from this fit. Adding a systematic 
error of 0.28 mag due to differential reddening to each data 
point would result P r (r 2 ) — 0.5 with 68% confidence ranges 
from 12.24 to 12.76 mag. We, therefore adopt a distance 
of 32 30+43?, pc for NGC 7419, which is similar to the es- 
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timate of ISubramaniam et alj (|2006|) but higher from the 
value report e d by Beauchamp. Moffat fc Drissenl (| 19941 ) and 
iBhatt et al.l l| 19931 ). 



4 NEAR-INFRARED DATA AND 
INTERMEDIATE MASS STARS 

To understand the global scenario of star formation around 
NGC 7419, we use NIR JHK 3 data for point sources within 
a square degree region centered on RA J20 oo = 22 h 52 m 28 s 
and DECj2ooo = +60°53'52", which is W.23 west from the 
cluster center. It contains the entire cluster and the sur- 
rounding star forming region. Data have been obtained from 
the 2MASS Point Source Catalog (PSC) (|Cutri et al.ll2003h . 
We have selected the sources based on the 'read-flag' which 
gives the uncertainties in the magnitudes. We retain 2MASS 
sources with 'read flag' values of 1-2 for good quality data 
in our analysis. The JHK S colors were transformed from 
2MASS system to California Institute of Technology system 
using the relations given on the 2MASS web site □ . 



4.1 Color-color diagrams 

The (J — H),(H — K) color-color diagram is plotted 
for the cluster region (471 stars) and the surrounding 
square degree region (22100 stars), in Figures [8] Solid 
lines represent unredde ned main-sequence and giant branch 
ijBessell fc Brettl Il988l ). The parallel solid lines are the 
reddening vectors for early and late-type stars (drawn 
from the base and t he tip of two branches). Location o f 
T-Tauri stars (TTS; iMever. Calvet fc Hillenbrand! 1 19971 ) . 
proto-star ( PS) like objects, class i cal Be stars and Her - 
big Ae/Be (|Doughertv et al.l Il994l : iHernandez et al.ll2005l ) 
are also shown. The extinction ratio Aj/Av = 0.282, 
Ah I Ay = 0.180 and Ak/Av = 0. 116 have been taken from 
ICardelli. Clayton fc Mathisl (|l989T l. Stars below the redden- 
ing vectors are considered to have NIR excess. The aster- 
isks represent all the stars with Ha emission (Table [6]). 
Of 44 emission-line stars only 39 stars have NIR coun- 
terparts within 3". NIR excess is shown by 17 emission- 
line stars within the cluster region but 4 more emission- 
line stars are situated at the boundary of the reddening 
vector, therefore they have also considered as NIR excess 
sta rs within error bars i n colo rs. From spectroscopic stud- 
ies, [Subramaniam~etjil] (2006) have identified them as Her- 
big Ae/Be stars. Therefore, NGC 7419 is a peculiar cluster 
having such a large number of Herbig Ae/Be stars. In sur- 
rounding region, 90 sources are fou nd in the T-Tauri locus 
ijMever. Calvet fc Hillenbrandlll997h . These sources are con- 
sidered to be mostly classical T-Tauri stars (Class II objects) 
with large NIR excesses and belong to the YSO population. 
There may be an overlap in NIR colors of Herbig Ae/Be 
stars and T-Tauri stars in the TTS region l|Hillenbrand et all 
1992). Such a large number of young stars represent the 
youth of this region. But all the T-Tauri stars with NIR 
excesses are situated outside the cluster region. 



http: / / www.astro.caltech.edu/~jmc/2mass /v3/ transformations / 
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Figure 9. Color-magnitude diagram My, (V — I)o of the cluster 
region. Open circles denote stars with Ha emission (Table O, 
asterisks denote stars with X-ray counterparts and having mem- 
bership probability from X-ray colors while filled dots represent 
stars with NIR excess. Post-main-sequence isochrone for 25 Myr 
by GRD02 (continuous line) and pre- main-sequence isochrones for 
0.3, 0.5, 1, 2, 4, 8 Myr by SES00 (dashed lines) are also shown. 
Isochrones are corrected for a distance of 3.23 kpc. 



4.2 Cluster Age from Herbig Ae/Be stars 

Figure shows My, (V — I)o color magnitude diagram for 
Ha emission stars (open circles; Table [6]), NIR excess stars 
(black dots), stars with X-ray counterpart (asterisks; H8.3p . 
Star BMD 950 is not shown as it shows a large error (> 
0.1 mag) in / band. It is seen that most of Herbig Ae/Be 
stars are located to the right side of the main-sequence and 
to es timate their age, we use pre-ma in-sequence isochrones 
from lSiess. Dufour fc Forestinil (|200fj| ). hereafter SES00, for 
ages 0.3, 0.5, 1, 2, 4 and 8 Myr, shown with dashed lines in 
Figure [9] For reference, the post-main-sequence isochrone of 
25 Myr age GRD02 have been shown by continuous lines. 
Only one X-ray source and three emission-line star lie in the 
range of 3.5-8 Myr isochrones, which may arise due to field 
star contamination of the Herbig Ae/Be stars loci in the NIR 
color-color diagram (Figure |SJ). Therefore, considering only 
the large number of Herbig Ae/Be stars, i.e. emission-line 
stars with NIR excess, it is very likely that the turn-on age 
of the cluster cannot be more than 2 Myr. 



4.3 Mass of Herbig Ae/Be stars using NIR 
color-magnitude diagram 

Stellar masses are determined using NIR color-magnitude 
J, (J — H) diagram. We prefer J over H or K, as the J 
waveban d is less affected by the emission from circumstellar 
material IjBertot. Basri fc Bouvierl 1988). Left panel of Fig- 
ure [10] shows NIR CMD for the cluster region NGC 7419. 
The 4 Myr post-main-sequence isochrone GRD02 and pre- 
main-sequence evolutionary tracks for mass range 2 — 7Mq 
SES00 are plotted assuming a distance of 3.2 kpc and a mean 
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Figure 8. Color-color diagrams using the 2MASS JHK data, left panel: The NIR color-color diagram for the cluster region NGC 7419. 
Small dots represent the NIR stars present in cluster region and the asterisks denote the Ha emission-line stars from Table [6] right panel 
: Same as the left panel, but for the square degree surrounding region centered at RA J20 oo = 22 h 52 m 28 s and DEC J2 ooo = +60°53'52". 
Bigger filled dots denote 90 T-Tauri stars in this region. 





(j-n) 
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Figure 10. Stellar mass estimate using NIR CMD diagram, left panel: J, (J — H) CMD for Herbig Ae/Be stars in the cluster region NGC 
7419 are shown with open circles. The 4 Myr post-main-sequence isochrone GRD02 and the pre-main-sequence evolutionary tracks for 
masses 2 — 7Mq are also shown and indicated, right panel: Same as left panel but for the YSO candidates (asterisks) in the surrounding 
square degree region of NGC 7419. Solid curve denotes pre-main-sequence isochrone of 6 Myr SES00, while the dashed oblique reddening 
lines denote the positions of pre-main-sequence stars for masses 0.1, 1, 2 and 3Mq. 



reddening E(B-V) of 1.70 mag ( ^3l^33]l . Reddening and 
extinction corrections to the isochrone and tracks are made 
using the relation Aj/Ay = 0.282 and A^j—u) /A v = 0.102 
l|Cardelli. Clayton fc Mathislll989l ), where A v = 3.1£(£ - 
V). Location of 21 Herbig Ae/Be stars in NIR CMD is rep- 
resented by open circles and it is seen that 66% of these are 
located between the mass range 5 — 7Mq . Only two stars lie 
in the mass range 3 — 4Mq . 



Right panel of Figure [10] represents NIR CMD for 90 
probable YSO candidates (shown with filled circles) iden- 
tified in a square degree region around the cluster (Fig- 
ure [8}. None of these YSOs are present in the cluster re- 
gion and hence we argue that the population of YSOs is 
most lik ely associated with a l ocal arm molecular cloud 
sh2-154 (|Ungerechts et al l I2000T ). see 35] for detailed dis- 
cussion. We, therefore, correct the theoretical isochrones 
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Table 9. IRAS point sources present in the square degree re- 
gion centered about RA J2 ooo = 22 h 53 m 26 s and DEC J2 ooo = 
+60°53'52'. ID 44 is a carbon star. 



ID 


IRAS PSC 


Fl2 


^25 


F 6 a 


F 100 


Remark 






(Jy) 


(Jy) 


(Jy) 


(Jy) 






22493+6018 


7.04 


2.39 


3.12 


48.91 


VMW Cop 


42 


22527+6030 


2.53 


0.77 


17.05 


87.48 


VMZCep 




22525+6033 


112.10 


92.30 


17.05 


50.10 


M7.5Iab 


14 


22520+6031 


3.10 


1.06 


17.05 


50.87 


M3.5 lab 




22512+6100 


108.50 


93.30 


16.68 


15.27 


V386Cep 




22538+6102 


1.04 


0.49 


6.50 


70.00 






22486+6028 


0.85 


0.37 


4.95 


53.07 






22490+6043 


0.54 


0.35 


5.35 


39.51 






22469+6053 


0.46 


0.25 


5.55 


71.76 






22480+6056 


0.30 


0.25 


5.34 


68.93 




Table 10. Same as 


for Table |9]except for the MSX point sources. 


ID 


Name 




F A 


F B 


F C 


Remark 




MSX PSC 


(Jy) 


(Jy) 


(Jy) 






G108. 6966+01. 0686 


6.70 


5.80 


3.72 


VMWCcp 




G109. 1599+01. 1206 


87.83 


133.58 


96.95 


M7.5 lab 


42 


G109. 1632+01.0645 


3.31 


2.14 


1.48 


VMZ Cep 


5 


G109. 1353+01. 1118 


3.40 


2.00 


1.68 


M2.5 lab 


4 


G109. 1302+01. 1321 


3.43 


2.11 


1.33 


M2.5 lab 


7 


G109. 1445+01.0947 


1.71 


0.71 


0.74 


M2.5 lab 


14 


G109. 0901+01. 1184 


3.38 


2.41 


1.58 


M3.5 lab 




G109. 2152+01. 6007 


63.48 


83.70 


66.06 


V386Cep 



SESOO, for extinction assuming the cloud distance of 1.4 kpc 
i|Khalil, Joncas fc Nekkall2004l ) and for reddening assuming 
a value of .4 mag/kpc towards the direction of NGC 7419 
l|joshill2005t ). The corrected isochrones are shown with solid 
curves in the right panel of Figure \W\ and they trace the lo- 
cus of 6 Myr pre-main-sequence stars having masses in the 
range 0.1 — 3Mq. It is observed that masses of all of the 
YSOs lie in the range 0.1 — 2.0Mq indicating them to be 
T-Tauri stars. Therefore, almost all of the YSO candidates 
are low mass pre-main-sequence stars. 



5 SPATIAL DISTRIBUTION OF LOW MASS 
PRE-MAIN-SEQUENCE CANDIDATES 

In Figure 1111 we show a square degree i?-band image, re- 
produced from DS^f|, of the field containing YSO candi- 
dates and the cluster NGC 7419. Inner (r < 1') and outer 
boundaries (r < 3.'5) of the cluster are encircled. Point 
sources from MSX (Midcourse Space Experiment) which 
surveyed the Galactic plane in four mid-infrared bands - 
A(8.28 /im), C(12.13 /im), D(1 4.65 /mi), and E( 21.34 pm) 
at spatial resolution of ~ 18" (|Price et al.ll2001j ). are also 
shown in Figure 1111 Information about 10 IRAS (Infra- 
Red Astronomical Satellite) point sources and 8 MSX point 
sources are given in Table [5] and Table 1101 respectively. 
The MSX sources are denoted by white open circles while 
asterisks represent IRAS sources. YSO candidates are de- 
noted by open square boxes. Of 8 MSX sources, six are 
identified with five supergiants and a carbon star and the 
other two are G108.6966+01.0686, a carbon star VMW Cep 
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Figure 11. A square degree i?-band DSS image of the field sur- 
rounding the cluster NGC 7419. North is up and east is to the left. 
Open boxes show the candidate YSOs, white thick dots represent 
MSX sources and open star symbols denote IRAS point sources. 
Inner (r < 1') and outer boundaries (r < 3.'5) of NGC 7419 
are encircled. The location of Sharpless nebula Sh2-154 is also 
marked. The dust extinction maps produced using star counting 
method are overlaid by white contours, and the ass ociated clumps 
P5, P 15 and P20 with cloud 699 are also marked IIDobashi et al.l 
120051) ■ 



ijKukarkin et al.lll968T ). and G 109. 2152+01.600 7, a semi reg- 
ular pulsating star V386 Cep (|Stephensonll984h . The cluster 
supergiants also have IRAS counterpart, and their presence 
in NGC 7419 is also supported from mid-infrared data. 

In the following, we describe the attenuation proper- 
ties of a square degr ee region around the cluster (see 
iDobashi et al.l l|2005h recently produced extinction maps of 
the entire region of the Galaxy in the galactic latitude range 
b < 40° by applying traditional star count technique to the 
optical DSS sources. We retrieved the FITS images of the 
extinction map of the field containing NGC 7419 and their 
surrounding region from their on-line websitfl In Figure [TT1 
we overlay the contours of high resolution extinction map. 
The contours are plotted at Av ~ 2.5 to 4.5 with an inter- 
val of 0.3 mag. The extinction towards the location of the 
cluster shows relatively low (Ay < 2.5) values, in compari- 
son to the outer region of the cluster especially, towards the 
west, the east and the south, where extinction increases up 
to ~ 4.5, ~ 3.2 and ~ 3.0 mag, respectively. These enhanced 
attenuations are identified with t he dense cloud clump s P5, 
P15 and P20 of dark cloud 699 l|Dobashi et al.l 120051 ). Us- 
ing star count method it would be difficult to detect dust 
obscuration located at a cluster distance of ~3.2 kpc, as 
the cloud would be come inconspicuous due to large numbe r 
of foreground stars ijDobashi et al.ll2005l : iMedhi et alJl2008h . 
We, therefore, argue that the dust obscuration would be due 
to a foreground cloud. 



Digital sky survey - http://stdatu.stsci.edu/dss/ 



http: / /darkclouds.u-gakugei. ac.jp/astronomer/astronomer. html 
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In order to see the dust emission characteristics of this 
cloud, we superimpose the IRAS 100 fim contours obtained 
from the Infrared Processing and Analysis Center (IPAC) on 
7?-band DSS image in Figure [T2l left panel. The IRAS survey 
was done in four bands 12 pim, 25 /im, 60 fim and 100 )im 
from mid-infrared to far-infrared. The contours are plotted 
after smoothing the IRAS 100 /im image. The contours are 
drawn at 134 (outermost) to 204 with the increment of 10 
M Jy sr~ . The IRAS dust emissions support the presence 
of dust obscuration at the location of cloud clumps P5, P15 
and P20. The star count extinction map is seen to be highly 
correlated with the 100 ^m IRAS dust map and represent 
the same morphology of the cloud. We also show the 12 CO 
temperature map of th e region, see right panel in Figure [121 
l|Kerton fc Bruntl2003l ) . The contours are drawn for temper- 
ature ranges 6.5 — 9.5 K with the increment of 0.5 K. The 
peak of CO emission is most likely associated with clump P5 
(~ 5' 5 westward). A slight mismatch might be the effect o f 
low resolution of the extinction map (|Dobashi et alj 12005). 
It is seen that 60% of total YSO population seems to be 
associated with the clumpy region P5. Therefore, the clump 
P5 may provide a fertile environment for the formation of 
the low mass stars. 

Interestingly, the spatial distribution of TT stars is 
highly correlated with the extinction and IRAS dust map 
(see Figure Hip , it is therefore likely that the population 
of YS Os (jjH) might be associated with the foreground dark 
cloud. lAvedisova (|2002l ) identify a local arm star forming re- 
gion Sh2-154, situated in the north-east direction from the 
cluster. Optical diameter and distance of Sh2-154 is esti- 
mate d respectively as 60' and 1.4 kpc l|Blitz. Fich fc Starkl 
Il982h . Hence, the YSO population (Figure [Sj) could be part 
of Sh2-154 and associated wi th a dark cloud clum p 699/P5 
and 699/P15 as identified bv lDobashi et all l|2005l) . 

The effect of environment on the production of the low 
mass stars has been investigated with the help of the distri- 
bution of the reddening in the cloud by assuming that they 
are situated at the same distance. The Ay value for each 
star was measured by tracing back to the intrinsic lines along 
the re ddening vector found in lMever. Calvet fc Hillenbrand! 
i| 19971 ). The value of Ay is found to vary from to 3 mag. 
Only one YSO has Ay ~5 mag. The spatial distribution of 
the Ay from to 5 mag is shown in Figure [131 left panel. No 
spatial gradient of Ay is found with the spatial distribution 
of the YSOs. The frequency distribution of Ay (Figure [TBI 
right panel) represents that 77% YSOs are having Ay less 
than 1.0 mag and 15% of YSOs with Ay in between 1.0 
to 2.0 mag. It shows that the environment for the YSOs is 
almost similar. 

The effect of the youthfulness of the YSOs in the spa- 
tial distribution can also be investigated with the help (H — 
^Oexcess, which is expected due to t he presence of th e cir- 
cumstellar disk dissipating with time lOasa et al. (120061). For 
YSOs , it is determined by the method used in Matsuvanagil 
(2006) and it lies in the range 0.002 - 0.3 mag. The dis- 
tribution of (H — /^excess (Figure 1141 left panel) does not 
show any spatial gradient with distribution of the YSOs, in- 
dicating that the YSOs having different age are distributed 
uniformly within the cloud. Frequency distribution of the 
H — K excess is shown in Figure 1141 right panel and it 
indicates that 89% of YSOs have nearly similar age. This 



represents a uniform distribution of the star formation rate 
within the cloud. 

The morphology of the cloud is almost the same for 
all the YSOs and we find very little age dispersion in the 
YSOs. There is no significant contribution of the external 
agents found in the formation of these YSOs in this cloud. 
Therefore, such a uniform distribution of YSO candidates 
might be the result of primordial fragmentation. 



6 MASS FUNCTION 

The initial mass function (IMF) is the distribution of stel- 
lar masses that form in a star formation event in a given 
volume of space. It is an important result of star forma- 
tion and together with the star formation rate, the IMF 
dictates the ev olution and fate of galaxies and star clusters 
!|Kroupall2002T ). However, the direct determination of IMF 
is not possible due to the dynamical evolution of stellar sys- 
tems. Therefore, we derived MF which is often expressed by 
the power law, N(logm) oc m r and the slope of the MF is 
given as : 



r = dlogN (logm) / dlogm 



(2) 



where N(logm) is the number of star per unit logarithmic 
mass inter val. In the solar neighborhood the classical value 
derived by ISalpeterl l|l955l ) is T = —1.35. In order to esti- 
mate the MF for the cluster region, we have used statistically 
cleaned sample described in £13.21 We have divided statisti- 
cally cleaned sample for two star formation episodes. The 
first star formation episode is at 25 Myr and the masses 
corresponding to V < 21 mag have reached to the main- 
sequence, therefore only main-sequence MF have been con- 
sidered for this episode. In Figure [TBI the MS represents the 
main-sequence band (selected in i]3.2[l and is drawn by long- 
dashed lines. The stars within this band are shown by solid 
dots. The main-sequence theoretical isochrone from GRD02 
is shown by solid line at 25 Myr age with upper limit at 
V ~ 15.0 mag because the stars above this point have been 
evolved from the main-sequence. Due to the above men- 
tioned reasons, we have considered the stars with V magni- 
tude ranges 15 — 21 mag for the main-sequence MF. The sec- 
ond episode is at age < 2 Myr, therefore pre-main-sequence 
MF has been derived for the stellar population attached with 
this episode. The pre-main-sequence theoretical isochrones 
by SES00 have been drawn by solid lines for ages 2 Myr 
and 0.1 Myr (upper limit of age in model), and by dashed 
lines with dots for masses ranges from 2 Mq to 7 Mq, re- 
spectively. The stars which are lying in between this age 
range have been used for the estimate of pre-main-sequence 
MF and shown by triangles in Figure [T51 The Herbig Ae/Be 
stars, described in |4] are shown by open circles. All these 
stars are excluded in the estimation of main-sequence MF 
even though they are lying inside the selected MS band, and 
considered as pre-main-sequence stars and included in the 
estimation of pre-main-sequence MF. We have used the stars 
shown by triangles in Figure [TJJ and the Herbig Ae/Be stars 
inside the MS band for the estimation of pre-main sequence 
MF for masses down to ~ 3Mq due to the incompleteness 
of mass below this limit. 

We have derived the MF from LF using the theoretical 
evolutionary models. The MF in two subregions i.e., the in- 
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Figure 12. Same as Figure ITT1 but white coun ters represent in left panel: Contours of 100 /im IRAS dust emissions, in right panel 
12 CO temperature map l lKerton fc Bruritll2003ri . 
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Figure 13. Distribution of the reddening, left panel : Spatial distribution of the reddening of the YSO candidates. Increasing sizes of 
filled boxes represent YSOs having Ay with ranges 0—1 mag, 1 — 2 mag, 2 — 3 mag and 4 — 5 mag, respectively. X-axis and Y-axis 
denote RAj2ooo an d DECj2000i respectively, right panel : Frequency distribution of the reddening for YSO candidates. 



ner (r < 1') and the outer (l'< r <3!5) region of the cluster, 
respectively, as well as for the whole cluster region (r < 3.'5) 
are given in Table [TT1 and plotted in Figure [TBI 

The slope of the MF of the main- sequence stars in the 
inner region of the clust er is estimated as T — —0.19 ± 0.15, 
which is natter than the lSalpeteil (|l955T ) value, while for the 
stars i n the outer reg ion, V = —1.52 ± 0.18, which is similar 
to the lSalpeterJ (|l955h value. The V value for the stars is thus 
steeper in the outer region indicating a preferential distri- 
bution of relatively massive stars towards the cluster center 
indicating a segregation of stellar masses. The slope for the 
MF for the whole cluster is V = —1.10 ± 0.19, which is sim- 
ilar to the Salpeter value in the mass range of 8.6 — 1AM q. 



This value of the MF is similar to t he value estimated by 
iBeauchamp. Moffat fc Drissenl (|l994l ) i.e. Y = -1.25 ± 0.10 
estimated at age~14 Myr and distance ~ 2.3 kpc for the 
cluster region. 



The slopes of the pre-main-sequence MF for the inter- 
mediate mass range 7 — 3A/0 in the cluster are derived as 
T = +1.54 ± 0.52 for inner, F = +1.64 ± 1.52 for outer, and 
T = +1.68 ± 0.84 for whole cluster region. The value of the 
pre-main-sequence MF slope is significantly flatter than the 
Salpeter value. 
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Figure 14. Age spread of YSOs. left panel : Spatial distribution of the (H — K) CXCCBB of the YSO candidates. The increasing sizes of 
filled boxes represent the YSOs having (H — K) CXCCBB in ranges — 0.05 mag, 0.05 — 0.10 mag, 0.10 — 0.15 mag, 0.15 — 0.20 mag and 
0.25 — 0.30 mag, respectively. X-axis and Y-axis denote RAj2ooo an d DECj2oo0i respectively, right panel: Frequency distribution of the 
(H — -ff)excess for YSO candidates. 



Table 11. The MF of the cluster NGC 7419. MS represent the MF for main-sequence population while PMS represents the MF for pre- 
main sequence population. The numbers of probable cluster members (N) have been obtained after subtracting the expected contribution 
of field stars in i]3.2l log<j!> represents log N(log m). 



Magnitude range 


Mass Range 


Mean 


Inner 


region 


Outer 


region 


Whole 


region 


(V mag) 


(Mq) 


log (M ) 


N 


log</> 


N 


log</> 


N 


logfli 


MS 


15.0 - 16.0 


8.59 - 7.10 


0.894 


9 


2.036 


6 


1.860 


15 


2.258 


16.0 - 17.0 


7.10 - 5.20 


0.789 


11 


1.912 


16 


2.075 


27 


2.302 


17.0 - 18.0 


5.20 - 3.57 


0.642 


14 


1.933 


27 


2.218 


41 


2.400 


18.0 - 19.0 


3.57 - 2.45 


0.479 


16 


1.991 


34 


2.318 


50 


2.486 


19.0 - 20.0 


2.45 - 1.78 


0.326 


13 


1.972 


65 


2.671 


78 


2.750 


20.0 - 21.0 


1.78 - 1.44 


0.207 


11 


2.172 


96 


3.008 


110 


3.067 


PMS 




7.00 - 6.00 


0.544 


4 


1.505 


7 


1.748 


11 


1.945 




6.00 - 5.00 


0.653 


3 


1.491 


11 


2.055 


11 


2.160 




5.00 - 4.00 


0.740 


5 


1.800 


4 


1.703 


9 


2.056 




4.00 - 3.00 


0.067 


5 


1.873 


16 


2.378 


21 


2.496 



7 MASS SEGREGATION 



There is evidence for mass segregation in a few Galactic 
as well as LMC star clusters, with the highest mass stars 
preferentially found towards the center of the cluster (see 
Saga r et al. Il988l; ISagar fc Richtlcr 1991; Sagar fc Griffiths 



Pandev et al.1 l200ll . 120051 : lKumar~Sa gar & Mclnick 



It is a well accepted finding that stars in the clus- 
ters evolve rapidly towards a state of energy equiparti- 
tion through stellar encounters i.e., mass segre gation. How- 
ever, observations of very young clus t ers (e.g., Sagar et all 
19881; iPandevJVIahra fc Sagarl Il992t iPandev et all 120051 ; 



Hillenbrand! 1 1997T ) suggest that the mass segregation may 
be the imprint of star formation itself. 



We have subdivided the main-sequence sample (used 
for the main-sequence MF in §5§ into two mass groups 
(3.6 < M/Mq < 8.6; 1.4 < M/M & < 3.6) to character- 
ize the degree of mass segregation in NGC 7419. Figure [T71 
shows cumulative distribution of main-sequence stars as a 
function of radius in two different mass groups. Effect of 
mass segregation can be seen in this figure in the sense that 
more massive stars (3.6 < M/M@ < 8.6) tend to lie towards 
the cluster center. The Kolmogorov-Smirnov (KS) test con- 
firms the statement that the cumulative distribution of mas- 
sive stars in the cluster is different from that of relatively 
less massive stars at a confidence level of 99%. There is a 
strong evidence of mass segregation in the main-sequence 
stars within this cluster. We estimated the relaxation time 
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Figure 15. Extinction corrected CMD at mean E(B — V)=1.70 
mag. Solid lines represent the theoretical isochrones for main- 
sequence at 25 Myr (GRD02) and pre-main-sequence at 2 Myr 
and 0.1 Myr (SES00). Solid dots within long-dash line show the 
stars selected for the estimation of main-sequence MF at 25 Myr 
i.e. first episode of star formation, while the triangles represent the 
stars used for pre-main-sequence MF. Herbig Ae/Be are shown 
by open circles. Herbig Ae/Be stars within main-sequence band 
(solid dots with open circles) have not been considered in esti- 
mation of main-sequence MF and included in the estimation of 
pre-main-sequence MF. 



to decide whether the mass segregation discussed above is 
primordial or due to dynamical relaxation. To estimate the 
dynamical relaxation time Te, we have used the relation 

8.9 x W 5 N 1/2 R h 3/2 



i 1 /2/ og (o.4iV) 



(3) 



where N is the number of cluster stars, Ru is the radius 
containing half of the clu ster mass and m is the average 
mass of the cluster stars (|Spitzer fc Hart|[l97ll ). The total 
number of stars in the cluster region are estimated as 321 
in the mass range (1.4 < M/Mq < 8.6). For the half mass 
radius, we have used half of the cluster extent i.e., 1.63 pc. 
Taking average mass of the cluster as 3.92 Mq, we have 
estimated the dynamical relaxation time Te for the cluster 
as 8.0 Myr, which is lower than the turn-off age of the cluster 
i.e., 19 - 25 Myr. Therefore, we can conclude that the cluster 
is dynamically relaxed. 
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Figure 16. MF determination of MF in the various subregions 
of the cluster NGC 7419. leftpanel : for main-sequence and, for 
pre-main-sequence. rightpanel : Inner region represents the core 
region (r < 1' ), Outer region represents the corona region (l'O < 
r < 3'5 ) and Whole region represents the cluster region (r < 3'5), 
as defined using radial density profile. The error bars represent 
1/y/N errors. 
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Figure 17. Cumulative radial distribution of main-sequence stars 
in mass intervals 1.44 < M/Mq < 3.57 and 3.57 < M/Mq < 
8.59. 



Table 12. XMM-Newton Observations of NGC 7419. 



8 X-RAY DATA ANALYSIS 

We have used the archival X-ray data from the XMM- 
Newton observations of NGC 7419. The observations were 
proposed by Christian Motch to search for Be+white dwarfs 
binaries in NGC 7419. The observation details are given in 
Table 1121 Our analysis is based on the CCD images from 
the European Photon Imaging Camera (EPIC). The de- 
tails of the X-ra y telescope and EPIC PN and MP S cam - 
eras are given by iJansen et al.l i|200lT ). IStriider et all (|200ll ) 
and lTurner et all ( 2001 ). Data reduction followed standard 



Parameter 



NGC 7419 



Observation ID 0201160501 
Start time (UT) 02 Feb 2004 02:40:39 

Stop time (UT) 02 Feb 2004 14:19:16 

Usable time(ks) 39.62(MOS1),40.21(MOS2),35.35(PN) 

EPIC mode full frame 

Optical filter Medium 
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Figure 18. The X-ray point sources in the mosaic image of MOS 
and PN detector in the energy band 0.3 - 7.5 keV with their ID 
number as given in Table [T3l X-axis and Y-axis denote RAj2ooo 
and DECj2000i respectively. 



procedures using the XMM-Newton Science Analysis Sys- 
tem software (SAS; version 7.0.0). Event files for MOS and 
PN are generated by using tasks emchain and epchain, re- 
spectively. Data from the three cameras were individually 
screened for high background episodes and the time inter- 
vals during which the total count rate (for single events of 
energy above 10 keV) in the instruments exceeded 0.35 and 
1.0 counts s _1 for the MOS and PN detectors, respectively, 
were excluded. 



8.1 Source detection and identification 

The source detection is based on the SAS point source de- 
tection algorithm edetectjzhain. Three energy ranges were 
selected - a soft (Sx) band (0.3 - 0.7 keV), a medium (Mx) 
band (0.7 - 1.2 keV) and a hard (Hx) band (1.2 - 7.5 keV), 
and we built the corresponding images for the different in- 
struments of the EPIC. Finally, source detection was per- 
formed on these images using edetectjzhain task, which is a 
chain script of various sub tasks. 

We inspected each source manually to reject false de- 
tections due to the instrumental artifacts, from the final list 
created by emldetect task. In this way, we found 66 sources 
with a combined maximum likelihood value greater than 10 
in all three instruments. The details of all the genuine X-ray 
sources in the energy band 0.3 - 7.5 keV are tabulated in 
Table [131 and shown in Figure [TBI 

We cross-correlated the X-ray source list with our opti- 
cal UBVRI phot ometric source list and NIR 2MASS All Sky 
data release data (|Cutri et al.ll2003D within the search radius 
6". The position of the stars in optical CCD are converted 
into RAj2ooo and DECj2ooo using the Guide Star Catalogue- 
II (CSC 2.2, 2001). There are 18 X-ray sources with opti- 
cal counterparts and 31 with NIR counterparts found after 



cross-correlation. Furthermore, none of the Herbig Ae/Be 
stars is having an X-ray counterpart. The unidentified X- 
ray sources i.e., without having any optical or NIR counter- 
parts within XMM-Newton field of view (30' x 30'), are 31 
in number (47% of the total) . 

8.2 The detection limit 

The analysis of detection limits in the XMM-Newton field 
of view is equally important to decide the X-ray emission 
level of undetected Herbig Ae/Be stars, as none of the Her- 
big Ae/Be stars is emitting X-rays. The faintest source de- 
tected above 2a in the XMM-Newton field of view (FOV) 
has a count rate of 8.4 x 10~ 4 counts s _1 in the PN detector. 
The detection limit is not constant throughout the field of 
view because of the following three reasons: (a) a gap be- 
tween the CCDs, and unoverlapped detector areas leading 
to a non-uniform effective detection area, (b) the decreasing 
exposure duration from the FOV center to its edges, causing 
a non-uniformity in the effective exposure times as well, and 
(c) the dense population of the stars near the center of the 
cluster affecting the detection limit non-uniformly. We have 
neglected the effects of the gap between CCDs, as an approx- 
imation, and investigated the effects of other parameters on 
the detection limit. 

We analyzed the exposure map created by the expmap 
task for the PN detector in the energy band 0.3 - 7.5 keV to 
derive the effect of exposure duration. It displays a smooth 
decrease in the counts from the center of the FOV to its 
edges by about a factor of three. If we are not considering 
the background effect, the signal-to-noise ratio is smaller by 
a factor of y/3 w 1.7 for the source at the edges of the FOV 
compare to the source at the center. Therefore, the detection 
limit in the edges of the FOV will be twice the detection 
limit at the center of the FOV due to the non-uniformity of 
exposure duration. 

The variations in the detection limit due to the dense 
clustering at the center of FOV ha ve been determined using 
the approach by I Sana et al.l l|2006l ). We computed an equiv- 
alent PN count rate in the energy band 0.3 - 7.5 keV due to 
the gap between CCDs. An empirical relation between count 
rates in the PN and MOS detectors was calculated and found 
to be approximately linear. For the sources which fall in the 
gaps between the PN CCDs, this relation is used to convert 
the MOS count rates into the PN count rates and these PN 
count rates are called pn-equivalent count rates. Figure [191 
left panel, displays the source pn-equivalent count rates as a 
function of the distance from the center of the FOV i.e., the 
center of the cluster NGC 7419. We adjusted a four-degree 
polynomial by selecting the faintest sources to derive an em- 
pirical detection limit in terms of pn-equivalent count rates 
(crum) as a function of distance (d) from the center of the 
FOV. This detection limit (in units of 10 -3 counts s _1 ) is 
shown by a solid line in Figure [191 left panel and described 
by the following relation: 

cr Um (d) = 2.246 - 0.999d + 0.267d 2 - 0.028d 3 + O.OOld 4 (4) 

where the distance d is expressed in arcmin from the cen- 
ter of the FOV. Left panel of Figure [19] shows that the 
detection limit is higher in center of the FOV (d < 3.'5), 
which is the cluster region as decided in £|3 . 1 1 It might be 
the effect of the extended X-ray emission from the cluster 
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Table 13. X-ray Sources detected using SAS task edetectjchain. Column 1 (XID), represents the identification number from X-ray 
source detection and display in Figure [18] Column 2 (Position), shows source position from X-ray source detection, Column 3 (Count 
Rates), represents the pn-equivalent count rates in energy band 0.3 - 7.5 keV, Column 4 (Hardness ratio), HR1 and HR2, defined in 
text, the sign of < represents the hardness ratio calculated after replacing zero counts by upper limits in cither of the energy band, 
Column 5 (Remarks), for PCM (the X-ray sources with membership probability using Hardness Ratio - the X-ray sources inside the 
dotted box in Figure I20H IR (the X-ray sources with NIR counterparts within 6"search radius using 2MASS data) , OPT (the X-ray 
sources with optical counterparts within 6" search radius using our UBVRI photometry.) , CI (the X-ray sources within cluster region 
NGC 7419 - see 



XID R--A-J2000 DECj2000 Count Rates Hardness ratio Remarks 

(h m s) (° ' ") counts s _1 HR1 H R2 



1 


22 


52 


35 


60 


47 


37 


10.92± 


1.50 


-0.40± 0.170 


0.34 ± 0.21 


IR 


2 


22 


52 


40 


60 


54 


06 


1.44± 0.58 


0.25 ± 0.393 


<-0.20 




3 


22 


52 


45 


60 


56 


01 


17.97± 


1.58 


0.16 ± 0.098 


-0.19 ± 0.11 


IR 


4 


22 


52 


46 


00 


54 


00 


2.61± 0.74 


-O.Olrb 0.276 


-0.58 ± 0.48 




5 


22 


52 


40 


60 


54 


57 


294. 66± 5.33 


0.47 ± 0.023 


-0.37 ± 0.02 


IR 


(i 


22 


52 


49 


60 


50 


01 


3.89± 


0.73 


0.82 ± 0.357 


0.28 ± 0.19 


IR, PCM 


7 


22 


53 


05 


00 


39 


19 


4.83± 


2.06 


<-0.43 


<0.64 




8 


22 


53 


00 


60 


57 


08 


4.15zt 


0.83 


0.19 ± 0.512 


0.77 ± 0.26 




9 


22 


53 


17 


00 


38 


13 


2.82± 


1.82 


0.54 ± 0.565 


<-0.63 


IR 


10 


22 


53 


21 


60 


41 


09 


5.05± 


0.73 


0.21 ± 0.167 


-0.10 ± 0.17 


IR 


11 


22 


53 


22 


60 


58 


08 


2.20± 


0.57 


0.43 ± 0.255 


<-0.43 




12 


22 


53 


20 


60 


40 


12 


3.91± 


0.54 


0.12 ± 0.145 


-0.34 ± 0.18 


IR 


13 


22 


53 


35 


60 


53 


12 


3.75:1 


0.59 


0.71 ± 0.683 


0.83 ± 0.21 


PCM 


14 


22 


53 


36 


60 


55 


47 


3.384: 


0.61 


0.69 ± 0.421 


0.54 ± 0.21 


PCM 


15 


22 


53 


37 


60 


44 


37 


0.844: 


0.31 


<0.64 


0.67 ± 0.44 


PCM 


16 


22 


53 


41 


61 


03 


43 


18.95± 


2.03 


0.40 ± 0.118 


-0.50 ± 0.14 


IR 


17 


22 


53 


42 


60 


59 


38 


3.18:1 


0.75 


0.04 ± 0.233 


-0.10 ± 0.31 


IR 


18 


22 


53 


52 


01 


05 


04 


6.12± 


1.14 


<0.97 


0.71 ± 0.23 


PCM 


19 


22 


53 


53 


60 


48 


48 


5.454: 


0.55 


0.83 ± 0.331 


0.70 ± 0.12 


IR, OPT, PCM, CI 


20 


22 


53 


57 


00 


59 


14 


1.32:1 


0.51 


0.58 ± 0.406 


-0.41 ± 0.46 


IR 


21 


22 


53 


57 


60 


56 


02 


8.794: 


0.96 


<0.85 


0.88 ± 0.15 


IR, PCM 


22 


22 


54 


00 


00 


54 


15 


1.564: 


0.41 


<0.70 


0.63 ± 0.31 


PCM 


23 


22 


54 


00 


00 


48 


22 


2.564: 


0.43 


<0.70 


0.75 ± 0.21 


OPT, PCM, CI 


24 


22 


54 


04 


60 


45 


45 


0.694: 


0.24 


0.08 ± 0.321 


0.13 ± 0.19 


IR, OPT 


25 


22 


54 


05 


60 


48 


04 


5.684: 


0.71 


0.07 ± 0.183 


0.48 ± 0.16 


IR, OPT, CI 


26 


22 


54 


07 


60 


43 


54 


1.754: 


0.37 


0.28 ± 0.230 


-0.22 ± 0.26 


IR, OPT 


27 


22 


54 


10 


60 


49 


29 


7.18:1 


1.17 


-0.56± 0.352 


0.84 ± 0.25 


CI 


28 


22 


54 


11 


60 


55 


23 


5.964: 


0.78 


0.30 ± 0.137 


-0.51 ± 0.18 


IR 


29 


22 


54 


12 


60 


47 


51 


1.644: 


0.37 


<-0.48 


0.75 ± 0.28 


IR, OPT, CI 


30 


22 


54 


12 


60 


45 


59 


1.624: 


0.33 


0.09 ± 0.20 


-0.98 ± 0.39 


IR, OPT, CI 


31 


22 


54 


10 


60 


50 


36 


5.614: 


0.73 


-0.16± 0.17 


0.41 ± 0.17 


IR, OPT, CI 


32 


22 


54 


18 


60 


48 


48 


4.78=1 


0.71 


<0.37 


0.53 ± 0.17 


IR, OPT, CI 


33 


22 


54 


21 


60 


59 


42 


24.48± 


1.52 


0.02 ± 0.06 


-0.65 ± 0.10 




34 


22 


54 


25 


60 


39 


07 


1.66± 


0.40 


0.49 ± 0.25 


<-0.25 


IR 


35 


22 


54 


27 


60 


48 


44 


1.40± 


0.38 


0.60 ± 0.52 


0.50 ± 0.31 


PCM, CI 


36 


22 


54 


29 


60 


43 


49 


1.64± 


0.55 


-0.06± 0.29 


<-0.18 


IR, OPT 


37 


22 


54 


29 


60 


35 


57 


7.26± 


2.37 


<-0.24 


<0.84 




38 


22 


54 


30 


60 


48 


12 


1.54± 


0.36 


0.85 ± 0.69 


0.65 ± 0.28 


OPT, PCM, CI 


39 


22 


54 


35 


60 


50 


05 


5.96± 


1.05 


0.60 ± 0.20 


-0.43 ± 0.21 


OPT, CI 


40 


22 


54 


39 


60 


41 


03 


4.48± 


0.59 


0.32 ± 0.14 


-0.58 ± 0.19 


IR. 


41 


22 


54 


42 


60 


44 


38 


3.33± 


0.84 


<-0.14 


<0.79 


IR, OPT 


42 


22 


54 


44 


60 


44 


01 


1.50± 


0.39 


<0.46 


<0.60 




43 


22 


54 


40 


60 


48 


27 


1.09± 


0.32 


0.37 ± 0.91 


0.80 ± 0.38 




44 


22 


54 


47 


60 


45 


18 


1.25± 


0.34 


<0.78 


<0.56 


PCM 


45 


22 


54 


49 


60 


52 


01 


1.14± 


0.35 


<0.78 


<0.53 


PCM 


46 


22 


54 


50 


60 


43 


29 


3.78± 


0.54 


0.93 ± 0.83 


0.87 ± 0.19 


OPT, PCM 


47 


22 


54 


55 


60 


47 


42 


2.25± 


0.42 


0.31 ± 0.21 


-0.03 ± 0.22 


IR, OPT 


48 


22 


54 


50 


60 


45 


37 


1.22± 


0.36 


0.43 ± 0.40 


0.23 ± 0.32 


IR, OPT 


49 


22 


54 


58 


00 


55 


27 


2.91± 


0.56 


<0.74 


0.88 ± 0.29 


PCM 


50 


22 


55 


02 


60 


53 


11 


1.40± 


0.48 


<0.86 


<0.60 


PCM 


51 


22 


55 


05 


60 


36 


45 


4.92± 


2.32 


<0.93 


<0.79 


PCM 


52 


22 


55 


10 


60 


54 


07 


3.07± 


0.60 


0.97 ± 0.72 


0.72 ± 0.24 


IR, OPT, PCM 


53 


22 


55 


17 


00 


40 


39 


0.92± 


0.37 


<0.70 


0.82 ± 0.52 


IR, OPT, PCM 


54 


22 


55 


22 


60 


53 


52 


3.12± 


0.61 


<0.96 


0.63 ± 0.23 


PCM 


55 


22 


55 


22 


60 


58 


10 


12.19± 


1.22 


0.01 ± 0.11 


-0.30 ± 0.14 


IR 


56 


22 


55 


24 


60 


47 


02 


3.35± 


0.97 


-0.23± 0.64 


0.76 ± 0.39 




57 


22 


55 


27 


60 


57 


18 


2.14± 


0.61 


0.20± 0.27 


-0.70 ± 0.44 




58 


22 


55 


30 


60 


41 


35 


4.14± 


1.05 


-0.62± 0.41 


0.74 ± 0.42 




59 


22 


55 


31 


60 


51 


39 


6.48± 


1.25 


-0.46± 0.45 


0.84 ± 0.28 




60 


22 


55 


34 


60 


42 


23 


1.85± 


0.55 


<0.99 


<0.54 


PCM 


61 


22 


55 


38 


60 


50 


42 


4.95± 


0.77 


<0.87 


<0.79 


PCM 


02 


22 


55 


44 


60 


50 


43 


5.48± 


1.08 


-0.76± 0.44 


0.91 ± 0.32 




63 


22 


55 


40 


60 


40 


48 


7.41± 


0.48 


<0.93 


<-0.12 


PCM 


04 


22 


55 


47 


60 


44 


58 


9.86± 


1.08 


0.58 ± 0.13 


-0.42 ± 0.13 


IR 


65 


22 


55 


54 


60 


47 


14 


2.22± 


0.63 


0.10 ± 1.50 


0.94 ± 0.39 




66 


22 


55 


50 


60 


56 


50 


2.94± 


0.83 


-0.02± 0.30 


0.14 ± 0.37 


IR 
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Figure 19. The variation of X-ray detection limits within XMM-Newton field of view in energy band 0.3 — 7.5 keV. leftpanel : The 
distribution of the pn-equivalent count rates of the X-ray sources as a function of their distance from the center of the cluster NGC 7419. 
The solid line shows the adopted lower limit given in equation (4). (b)rightpanel: Lower panel : estimated detection limit expressed 



in terms of the observed flux (log(erg s 
distance of the cluster NGC 7419. 



2 )), Upper panel : the detection limit, expressed in luminosity (log(erg s 1 )) using the 



region, discovered in §83] The detection limit increases to- 
wards the CCD edges, as estimated by the exposure map. 
Using mission count rate simulator, WebPIMMifl we con- 
verted the pn-equivalent count rates into fluxes and lu- 
minosities. Assuming a hydrogen column density (Nh) of 
9.9 x 10 21 cm" 2 , at E(B - V) =1.70 mag for the cluster 
NGC 7419, the Raymond-Smith single temperature opti- 
cally thin thermal plasma models were used at tempera- 
tures (kT) of 0.2 keV, 1.0 keV and 2.0 keV. Results are 
displayed in Figure 1191 right panel for different plasma tem- 
peratures. Luminosities are calculated for the distance of 
the cluster NGC 7419 i.e., 3.2 kpc. Within the cluster region 
(d < 3'5) , the detection limit of flux varies from 1.66 x 10 
to 3.89 x 10" 15 erg cm" 2 s _1 for soft X-ray sources at kT=0.2 
keV, 3.02 x 10" 15 to 6.31 x 10" 15 erg cm -2 s _1 for X-ray 
sources at kT=1.0 keV, and 4.47 x 10~ 15 to 10.0 x 10" 15 
erg cm -2 s" 1 for X-ray sources at kT=2.0 keV. We have 
estimated the detection limit of 21 Herbig Ae/Be stars, 
according to their positions in the cluster region using 
equation 3] As majority of Herbig Ae /Be stars have X- 
ray temperatures greater than 1.0 keV (jStelzer et al.l 120061 ; 
lHamaguchi. Yamauchi fc Kovamall2005r ). the detection lim- 
its for these stars are estimated as the median of the limiting 
X-ray luminosity (Lx) of Herbig Ae/Be stars i.e., 5.20 x 10 30 
erg s _1 . 



8.3 Probable Members Using Hardness Ratio 

We have determined the X-ray hardness ratio HRl and HR2 
defined as 



7 WebPIMMS is a NASA's HEASARC tool 
ered by PIMMS version 3.9. It is hosted at 
http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



pow- 
URL: 



HRl 



HR2 = 



(Mx+Sx) 



[H> 



My 



(Hx + M x ) 



(•») 



(6) 



where S x , Mx and Hx denote the soft band (0.3 - 0.7 keV), 
the medium band (0.7 - 1.2 keV) and the hard band (1.2 - 
7.5 keV), respectively. For undetected sources, we used up- 
per limits as a function of distance from the center of FOV 
to its edges in each the energy bands using the method men- 
tioned in £18.21 In cases where no counts are observed in any 
one energy band, the HRl or HR2 axe either +1.0 (no soft 
counts) or -1.0 (no hard counts), we have replaced the zero 
counts value by the upper limits depending upon the source 
position in the CCD in that energy band. The hardness 
ratios HRl and HR2 of the X-ray sources in the XMM- 
Newton FOV are displayed in Figure [201 X-ray sources are 
represented by solid dots, X-ray sources having optical coun- 
terparts, NIR counterparts and extended behavior are de- 
noted by open circles, open squares and the symbol of star, 
respectively. X-ray sources with rightward arrows represent 
the upper limits in HRl, and downward arrows represent 
the upper limits in HR2. Solid lines represent hardness ra- 
tios derived from model spectra. 

We have simulated the values of HRl and HR2 us- 
ing plasma model APEC (Atomic Plasma Emission Code) 
to relate these hardness ratios to the spectral properties 
of the sources and to decide the membership of X-ray 
sources in the cluster region. Using the EPIC-pn response 
matrices within XSPEC, we have generated model spectra 
for monothermal plasma (IT- APEC) and two-temperature 
plasma (2T-APEC : the coronal X-ray emission from active 
late-type stars is generally not monother mal and consistent 
with the 2T thermal p l asma models, (see | Favata et ai1l2003l ; 
iTsuiimoto et all |2002| ; IStelzer et all |2006| )). We considered 
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Figure 20. HR1 vs HR2 diagram for the X-ray sources denoted by solid dots. X-ray sources with optical counterparts, NIR, counterparts 
and extended behavior are denoted by open circles, open squares and the symbol of star, respectively. X-ray sources with rightward arrow 
represent the upper limits in HRX, and downward arrow represent the upper limits in HR2. Solid lines represents HR1 vs HR2 curves 
obtained from APEC model and dotted box contains the X-ray probable members sharing the cluster environment, leftpanel For the 
single temperature 1T-APEC plasma model, we considered plasma temperatures of kT ranges from 0.2 to 8.0 keV, at different hydrogen 
column density (N H )= 1.0 X 10 20 , 1.0 X 10 21 , 4.9 X 10 21 , 9.9 X 10 21 (for cluster region), 19.8 X 10 21 cm -2 . The temperature increases 
from bottom to top for each Nu value, rightpanel The two temperature 2T-APEC plasma model with kTl = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8 keV with the combination of kT2 ranges from 0.9 to 9.6 keV, at different hydrogen column density (Nh)= 1-0 X 10 21 , 4.9 X 10 21 , 
9.9 X 10 21 (for cluster region), 19.8 X 10 21 cm . The temperature kTl increases from left to right for each Njj value with increasing 
kT2 from bottom to top. 



plasma temperatures (kT) from 0.2 keV to 8.0 keV for 1T- 
APEC model, and kTl from 0.2 to 0.8 keV with the combi- 
nation of kT2 from 0.9 to 9.6 keV for 2T-APEC model at dif- 
ferent hydrogen column density (A?n)=1.0 x 10 20 , 1.0 x 10 21 , 
4.9 x 10 21 , 9.9 x 10 21 (for cluster region), 19.8 x 10 21 cm" 2 . 

The X-ray sources lying in between 1 x 10 20 < Nh < 
1 x 10 21 cm -2 are probably the foreground stars having ei- 
ther NIR counterparts or optical counterparts. As the value 
of mean Nh for the cluster is estimated as 9.9 x 10 21 cm" 2 
from the optical studies, X-ray stars located in between 
4.9 x 10 21 < N H < 19.6 x 10 21 cm" 2 belong to nearly the 
same cluster environment. In Figure 1201 the X-ray sources 
situated inside the dotted box are considered as probable 
members belonging to the cluster. Most of the X-ray sources 
inside this box are having only upper limits in the soft bands 
and have hard spectra, which is a consequence of the ab- 
sorption of the soft components of the energy by the large 
column densities. A few sources are fitted neither by 1T- 
APEC model nor 2T-APEC model and show a soft energy 
component as well as a hard energy component. Most of 
these sources are unidentified sources not having any NIR 
or optical counterparts. They may be the foreground stars, 
whose positions are uncertain in hardness ratio diagram be- 
cause of the large uncertainty in the hardness ratios. But, 
for a few X-ray sources, the uncertainty in hardness ratios is 
not so large that their location can be explained within the 
temperature range of the models. They may be either the 
foreground stars having very high hard energy component 
greater than 9 keV or the IT- APEC and 2T-APEC model 
are not enough to characterise them. 
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Figure 21. XLF for the probable members in the clus- 
ter region NGC 7419, weak-lined (WTTS) and classical 
T-Tauri stars (CTTS) in Taurus- Auriga regio n, and Her- 
big Ae/Be (HAe/Be) stars from | Stelzer et all l|2006h and 
lHamaeuchi, Yam auchi & Koyama ( 20Q3) using Kaplan-Meier es- 
timator. The downward arrows represent 21 Herbig Ae/Be stars 
in the cluster NGC 7419. 



8.4 X-ray Luminosity Function 

X-ray luminosity function (XLF) is frequently employed 
to characterize a stellar population. We performed statis- 
tical analysis of XLF using Kaplan-Meier estimator of in- 
tegral distribution functions. We have derived the median 
Lx of weak-lined and classical T-Tauri stars in Taurus- 
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Figure 22. The X-ray extended emission from the cluster NGC 
7419 in mosaic MOS detector in energy band the 1.2 - 7.5 keV. 
X-axis and Y-axis denote RAj2ooo and DECj2000i respectively. 
The gray scale is logarithmic to highlight the diffuse component. 
The contours are plotted at 4.104 (above 3<r of the background), 
4.644, 5.418, 6.192, 6.966 10 -3 counts s _1 arcrnin" 2 . 



Auriga region (age= 0.1 — 10 Myrs) using the data from 
IStelzer fc Neuhauserl l|200l[ ) and for Herbig A e/Be stars, 
we have used the sample of Stelzer et al.l l|2006h and 
lHamaguchi. Yamauchi fc Kovamal 1 2005t ). The median Lx 
of four probable members in the cluster is estimated as 
30.75 erg s -1 . For weak-lined T-Tauri, classical T-Tauri, and 
Her big Ae/Be stars in the sample of Istelzer et al.l (|2006l ) 
and lHamaguchi. Yamauchi fc Kovamal (|2005l ). the median 
L x have been estimated as 29.88, 29.07, 29.29 and 31.15 erg 
s -1 , respectively. A comparison of the XLFs of T-Tauri stars 
and Herbig Ae/Be stars with the cluster members are shown 
in Figure [2T1 

The X-ray upper limits indicate that the Herbig Ae/Be 
stars in NGC 7419 are not systematically more active than 
T-Tauri in Taurus- Auriga region and Herbig Ae/Be stars in 
IStelzer et all (I2006T). But the median value of XLF of Her big 



Ae/Be stars in Hamaguchi. Yamauchi fc Kovamal l|2005t ) is 
more than the upper limits in the cluster. Our exposure and 
resolution (see t|8.5jl are not enough to reach the median 
value of Lx of T-Tauri stars, therefore, we can not conclude 
whether Herbig Ae/Be stars can have an XLF similar to ei- 
ther t he T-Tauri stars or Herbig Ae/Be stars in lStelzer et alj 
(2006). A much deeper detection level is clearly needed in 
order to reach definitive conclusions concerning in the X- 
ray emission from Herbig Ae/Be stars. However, the more 
important conclusion to draw from this analysis is that the 
Herbig A e/Be stars in NGC 7419 are not as a ctive as in the 
sample of lHamaguchi. Yamauchi fc Kovamal (|2005T ). 



8.5 Extended X-ray Emission 

The extended emission from the cluster region NGC 7419 
has been studied using MOS images. PN data were not used 
since chip gaps, bad columns, and the the out-of-time event 
correction would have complicated the analysis. We have 
created the MOS images and associated exposure maps in 
three bands: soft band (0.3 - 0.7 keV), medium band (0.7 




Figure 23. Overlay of X-ray point sources (black open circles) 
with ID from Table IT3l and 21 Herbig Ae/Be stars (white open 
circles) in the cluster region NGC 7419. Large black circle repre- 
sent the cluster region and black squares represent supergiants. 
X-axis and Y-axis denote RAj2ooo and DEC.J200O; respectively. 
The white contours represent the extended X-ray emission in en- 
ergy band 1.2 - 7.5 keV above the 3<r of the background count 
rate i.e. the outermost contour of the Figure l22l 



- 1.2 keV) and hard band (1.2-7.5 keV). We have cho- 
sen adaptive smoothstyle with signal-to-noise ratio 10 in the 
task asmooth, for getting the faint emission from the clus- 
ter region. The resulting mosaic image with point sources 
are weighted, masked and exposure-corrected following the 
procedure described in the documentation for the SAS task 
asmooth^. The point sources are removed from the image 
using cheesemask but the background is not subtracted from 
it. 

The signature of extended emission is not present in 
the soft band (0.3 - 0.7 keV) and medium band (0.7 - 
1.2 keV) images of the MOS, because of the the absorp- 
tion of the soft X-ray component due to the large hydro- 
gen column density (Nh)- All the probable point sources 
have hard spectra, therefore, it is very likely that the ex- 
tended emission is due to unresolved hard X-ray sources 
in the cluster. The X-ray contour map of the diffuse X-ray 
emission in energy band 1.2-7.5 keV is shown in Figure |2"21 
The entire cluster region is covered by the central CCD of 
the MOS, therefore, we have chosen background regions of 
22745 pixels beyond the diffuse emission but within the cen- 
tral CCD. The contours are plotted above the 3<r of the 
mean background count rates. The extended X-ray emission 
from the cluster region contains 0.084 counts s -1 in 14611 
pixels and corresponds to 5. 76x 10 -6 counts s -1 pixel -1 . 
The pixel size is set at 2" x 2". Therefore, the total ex- 
tended emission has a size of 16.23 arcmin 2 . We have 
estimated the count rate from the extended emission as 
1.38xl0 -3 counts s -1 arcmin -2 , after subtracting the mean 
background count rate. Using WebPIMMS ver3.9a, this 
count rate is converted into flux after considering Raymond- 
Smith thermal plasma model at 1 keV and estimated as 



1.47 x 10 



-14 -1 

erg s cm 



arcmin , which corresponds to 



http: / /xmm. gsfc.nasa.gov/docs/xmm/sas/help/asmooth/index. html 
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the luminosity Lx ~ 1.8 x 10 31 erg s _1 arcmin -2 in the en- 
ergy band 1.2 - 7.5 keV. We have chosen the beamwidth as 
50 pixel i.e., 5er (smoothing beam width). We have estimated 
18 beam elements in per arcmin 2 region and estimated Lx 
of unresolved sources as 1.0 x 10 30 erg s" 1 . Therefore, wo 
are expecting ~ 288 unresolved X-ray sources in the total 
diffuse component. 

The Lx of unresolved X-ray sources is comparable with 
the Lx of T-Tauri stars i.e., 28 < log Lx < 32 in energy 
band 0.5 to 8.0 keV , with a peak around log Lx ~ 29, 
|Feigelson et alj|2005h . It can be explained by the presence of 
~ 288 T-Tauri stars in the cluster, which are not resolved by 
XMM-Newton observatory. Higher resolution (sub arcsec) 
data are urgently required to explore this issue. 



8.6 X-Ray emission and Herbig Ae/Be stars 

We have examined a sample of 21 Herbig Ae/Be in this 
cluster using 2MASS data, H a photometry and the spec - 
troscopic studies reported by ISubramaniam et all (2006). 
However, none of Herbig Ae/Be stars has an X-ray coun- 
terpart. In Figure 1231 the distribution of the X-ray point 
sources which are having optical counterparts, supergiants 
(5 in number) and the Herbig Ae/Be stars within the clus- 
ter are shown by black open circles (with ID from Table [T3)). 
black open squares and white open circles, respectively. 

The detection limit analysis £18.21 and the XLF (Fig- 
ure I21|) indicate that the Lx observed from the low mass 
pre-main-seq uence stars is less than upp er limits of the de- 
tection (see, IStelzer fc Neuhauserl 120011 ) i.e., 28 < log L x 
< 32 in energy band 0.5 to 8.0 keV, with a peak around 
log Lx ~ 29 erg s" 1 , (Feigelson et al. 2005). Therefore, if a 
Herbig star is having a T-Tauri star as a binary companion, 
then we cannot detect the X-ray emission from it. There- 
fore, the generation of X-ray emission from Herbig Ae/Be 
stars (see, Zinnecker fc Preib isch 1994:; ll5amiani et al.ll 1994 ; 
IStelzer et alj|2006h might be the result of a T-Tauri binary 
companion which is not detectable in our study and we can 
not rule out the companion hypothesis for the generation of 
X-ray emission from Herbig Ae/Be stars. It might happen 
that the Herbig Ae/Be star itself is emitting X-rays but the 
level of the X-ray emission is less than the detection limit. 
A similar kind of process as in T-Tauri stars could then take 
place in the Herbig Ae/Be stars for the generation of X-rays. 



9 SUMMARY AND CONCLUSIONS 

A deep optical UBVRI and narrow band Ha observations 
along with multi-wavelength archival data from the surveys 
such as 2MASS, MSX, IRAS and XMM-Newton are used 
to understand the global scenario of star formation and the 
basic parameters of the cluster NGC 7419. XMM-Newton 
archival data have also been used to study the X-ray emis- 
sion mechanism from the cluster. 

The radius of the cluster NGC 7419 has been found 
to be 4f0±0f5 using radial density profile. The reddening 
law in the direction of the cluster is found to be normal at 
longer wavelengths but anomalous at shorter wavelengths. 
Reddening, E(B — V), is found to be varying between 1.5 
to 1.9 mag with a mean value ~ 1.7 ± 0.2 mag. The turn- 
off age and the distance of the cluster are estimated to be 



22.5 ± 2.5 Myr and 3230±|^, respectively. The MF for the 
main-sequence stars in the cluster is estimated as having 
r = -1.10 ± 0.19 in th e mass range 8. 6M < M < 1.4M , 
which is a similar to the lSalpeterl (|1_955) value. Effect of mass 
segregation is found in the main-sequence stars which may 
be the result of dynamical evolution. 

Using the NIR color-color diagram and narrow band 
Ha observations, we have identified 21 Herbig Ae/Be in the 
cluster region with the masses lying between 3 to 7Mq. The 
ages of these Herbig Ae/Be stars are found to be in the range 
of ~ 0.3 to 2.0 Myr. The significant difference between turn- 
off age and turn-on age of the cluster represents a second 
episode of star formation in the cluster. We have found 90 
YSOs having masses in the range from 0.1 to 2.OM0 with 
the help of NIR color-color diagram around the cluster. The 
presence of such a large number of NIR excess sources (T- 
Tauri stars) shows a recent star formation episode in the 
surroundings of the cluster region. Using extinction, dust 
and 12 CO maps, we found that these YSOs are probably 
associated with a foreground star forming region Sh-154 and 
not related with the cluster region. We found no obvious 
trend in spatial distribution of YSOs with Ay and (H — 
^Oexcess- The dispersion in Ay and (H — K) axccas is also 
very low which indicates that a majority of the YSOs are 
born at the same time in the same environment. Therefore, 
it is possible that primordial fragmentation of the cloud may 
be responsible for the formation of the low mass stars within 
this cloud. 

We have detected 66 X-ray sources in the observed field 
by XMM-Newton observatory using archival X-ray data. 
Out of these sources, 23 are knowm to be the probable mem- 
bers of the cluster based on analysis of their X-ray colors. 
Fifteen X-ray sources are without any optical or NIR coun- 
terparts. These may be the young embedded sources which 
need to be investigated further. We have derived the detec- 
tion limits of X-ray observations based on the position of the 
X-ray source in the field of view and on its energy spectrum. 
Thus, the median value for the detection limit for the 21 
Hebig Ae/Be stars in the field is L x ~5.2 x 10 30 erg s -1 . We 
have compared the XLF for the cluster members with the T- 
Tauri stars and Herbig Ae/Be stars. Because of insufficient 
exposure and resolution, the sensitivity was not enough to 
reach the level of the median X-ray luminosity observed in T- 
T auri stars in Taurus- A uriga region and Herbig Ae/Be stars 
in IStelzer et al.l (|2006T ). Therefore a conclusive comparison 
of X-ray properties of the stars cannot be made. However, 
the comparison indicates that Herbig Ae/Be stars in NGC 
7419 tend to be less X-ray luminous tha n in the sample of 
lHamaguchi. Yamauchi fc Kov ama ( 2005) , which shows that 
X-ray activity level of the Herbig Ae/Be stars is not more 
than in the T-Tauri stars. Therefore, we can support the bi- 
nary T-Tauri companion hypothesis for the generation of X- 
rays in Herbig Ae/Be stars. It is also possible that a Herbig 
Ae/Be star is itself emitting X-rays but the level of the X-ray 
emission is similar to that of the T-Tauri stars. The cluster 
region shows an extended X-ray emission with a total lu- 
minosity estimated to be Lx w 1.8 x 10 31 erg s _1 arcmin -2 . 
This diffuse emission might be the result of X-ray emis- 
sion from T-Tauri type stars which could not be resolved by 
XMM-Newton observations. It requires ~288 T-Tauri stars 
each having Lx ~ 1.0 x 10 30 erg s _1 , if it originates from 
such stars. High resolution deep observations such as from 
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CHANDRA are required for a detailed analysis of this clus- 
ter region. 
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